Technical  Report  EL-98-5 
June  1998 


US  Army  Corps 
of  Engineers 

Waterways  Experiment 
Station 


Use  of  Nutrient  Balances  in  Comprehensive 
Watershed  Water  Quality  Modeling 
of  Chesapeake  Bay 

by  Anthony  S.  Donigian,  Jr.,  Radha  V.  Chinnaswamy 
Aqua  Terra  Consultants 

Patrick  N.  Deliman,  WES 


Approved  For  Public  Release;  Distribution  Is  Unlimited 


1 99807 OT  225 


DTIC  QUAUTy  mSPECTED  1 

Prepared  for  U.S.  Environmental  Protection  Agency 


The  contents  of  this  report  arc  not  to  be  used  for  advertising, 
publication,  or  promotional  purposes.  Citation  of  trade  names 
docs  not  constitute  an  official  endorsement  or  approval  of  the  use 
of  such  commercial  products. 

The  findings  of  this  report  arc  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 


PRINTED  ON  RECYCLED  PAPER 


Technical  Report  EL-98-5 
June  1 998 


Use  of  Nutrient  Balances  in  Comprehensive 
Watershed  Water  Quality  Modeling 
of  Chesapeake  Bay 


by  Anthony  S.  Donigian,  Jr.,  Radha  V.  Chinnaswamy 

Aqua  Terra  Consultants 

2672  Bayshore  Parkway,  Suite  1001 

Mountain  View,  CA  94043-1001 

Patrick  N.  Deliman 

U.S.  Army  Corps  of  Engineers 
Watenvays  Experiment  Station 
3909  Halls  Ferry  Road 
Vicksburg,  MS  39180-6199 


Rnal  report 

Approved  for  public  release;  distribution  is  unlimited 


Prepared  for  Chesapeake  Bay  Program  Office 

U.S.  Environmental  Protection  Agency 
Annapoiis,  MD  21403 


US  Army  Corps 

of  Engineers 

Waterways  Experiment 
Station 


NEAOOUWnERS 


ENV»KMMP(TM. 
MOORATOnY  — 


FOR  INFORMATION  CONTACT; 

PUBLIC  AFFAIRS  OFFICE 
U  S.  ARMY  ENGINEER 
WATERWAYS  EXPERIMENT  STATION 
3909  HALLS  FERRY  ROAD 
VICKSBURG,  MISSISSIPR  39180-6199 
PHONE:  (601  >634-2502 


STRUCTURES 

UBOfMTORr 


AREA  OF  RESERVATKW  .  2  r  * 


Waterways  Experiment  Station  Cataloging-in-Publication  Data 
Donigian,  Anthony  S. 

Use  of  nutrient  balances  in  comprehensive  watershed  water  quality  modeling  of 
Chesapeake  Bay  /  by  Anthony  S.  Donigian,  Jr.,  Radha  V.  Chinnaswamy,  Patrick  N.  Deliman 
;  prepared  for  Chesapeake  Bay  Program  Office,  U.S.  Environmental  Protection  Agency. 

1 18  p. :  ill. :  28  cm.  -  (Technical  report ;  EL-98-5) 

Includes  bibliographic  references. 

1.  Water  quality  -  Mathematical  models.  2.  Environmental  quality  -Mathematical 
models.  3.  Valley  ecology  -  Mathematical  models.  4.  Stream  ecology  -  Mathematical 
models.  4.  Chesapeake  Bay  (Md.  and  Va.)  -  Mathematical  models.  I.  Chinnaswamy, 
Radhakrishnan  V.  II.  Deliman,  Patrick  N.  III.  United  States.  Army.  Corps  of  Engineers.  IV. 
U.S.  Army  Engineer  Watenways  Experiment  Station.  V.  Environmental  Laboratory  (U.S. 
Army  Engineer  Waterways  Experiment  Station)  VI.  United  States.  Environmental  Protection 
Agency.  Chesapeake  Bay  Program.  VII.  Title.  VIII.  Series:  Technical  report  (U.S.  Army 
Engineer  Waterways  Experiment  Station) ;  EL-98-5. 

TA7  W34  no.EL-98-5 


Contents 


Preface  . vi 

Conversion  Factors,  Non-SI  to  SI  Units  of  Measurement .  vii 

1- Introduction .  1 

Study  Background .  1 

Scope  and  Objectives . 1 

Summary  Conclusions  and  Recommendations  .  2 

Format  of  this  Report  .  6 

2- Expected  Nutrient  Balances  by  Land  Use  Categoric .  7 

Overview  and  Summary  Tables .  7 

Forest  Balances .  8 

Pasture  Balances .  8 

Urban  Balances .  9 

Agricultural  Cropland  Balances .  11 

Qosure .  1 1 

3- Application  ofNutiient  Balances  and  Model  Refinements  to  the  Shenandoah  River  Watershed . 14 

Overview  of  Model  Refinements  .  14 

Hydrology  Simulation .  15 

Nonpoint  Source  and  LoacSng  Assessment .  16 

Water  Quality  Calibration  Results .  25 

Conclusions  and  Recommendations .  28 

References .  38 

Appendix  A:  Shenandoah  Model  Segments  Results .  A1 

SF  298 


iii 


List  of  Figures 


Figure  1. 
Figure  2. 

Figure  3. 
Figure  4. 
Figure  5. 
Figure  6. 

Figure  7. 

Figure  8. 

Figure  9. 

Figure  10. 

Figure  11. 

Figure  12. 

Figure  13. 


Chesapeake  Bay  Watershed  Model  and  Shenandoah  Model  Segments  .  3 

Simulated  and  Observed  Flow  at  Reach  190,  SF  Shenandoah  River 

at  Front  Royal,  VA .  18 

Simulated  and  Observed  Flow  at  Reach  200,  Shenandoah  River  at  Millville,  WV .  19 

Frequency  Analysis  of  Flow  at  Reach  190,  SF  Shenandoah  River  at  Front  Royal,  VA  ...  20 

Frequency  Analysis  of  Flow  at  Reach  200,  Shenandoah  River  at  Millville,  WV .  21 

Simulated  and  Observed  Sediment  (TSS)  Concentration  at  Reach  200,  Shenandoah  River 
at  Millville,  WV .  30 

Simulated  and  Observed  NO3-N  Concentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  31 

Simulated  and  Observed  NH3-N  CcHicentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  32 

Simulated  and  Observed  Organic  N  Concentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  33 

Simulated  and  Observed  Total  N  Ccmcentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  34 

Simulated  and  Observed  PO4-P  Concentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  35 

Simulated  and  Observed  Organic  P  Concentration  at  Reach  200,  Shenandoah  River  at 
Mmvi]le,WV .  36 

Simulated  and  Observed  Total  P  Concentration  at  Reach  200,  Shenandoah  River  at 
Millville,  WV .  37 


IV 


List  of  Tables 


Table  1.  Typical  Nitrogen  Balances  for  Major  Crops  and  Land  Use/Land  Cover  Categories  . 12 

Table  2.  Typical  Phosphorus  Balances  for  Major  Crops  and  Land  Use/Land  Cover  Categories _ 13 

Table  3.  Shenandoah  River  Watershed  Hydrologic  Calibration:  Comparison  of  Annual  Total 

Observed  versus  Simulated  Flow . .  17 

Table  4.  Unit  Area  Norqx)int  Source  Loading  Rates  for  Each  Land  Use 

for  the  Shenandoah  Basin .  23 

Table  5.  Percent  of  Total  Load  Contributed  Frcsn  Each  Source  in  Shenandoah  Basin .  23 


V 


Preface 


The  work  herein  was  authorized  under  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
Contract  No.  DACW39-94-C-0052  with  Aqua  Terra  Consultants,  dated  30  March  1994  and  amended 
25  July  1995;  the  work  was  completed  on  30  September  1996. 

Dr.  Patrick  N.  Deliman,  Environmental  Laboratory  (EL),  WES,  was  the  Project  Officer  and  co-author 
of  this  report.  He  was  instrumental  in  facilitating  the  administration  and  execution  of  the  contract  work. 

The  Use  of  Nutrient  Balances  in  Comprehensive  Watershed  Water  (Quality  Modeling  of  Chesapeake 
Bay  study,  as  documented  in  this  report,  was  performed  for  the  U.S.  Environmental  Protection  Agency 
Chesapeake  Bay  Program  Office  (CBPO),  Aimapolis,  MD.  Mr.  Lewis  Linker,  CBPO,  was  point  of 
contact  The  CBPO  provided  data  for  model  testing  and  refinement  within  the  Chesapeake  Bay 
Watershed  Model  that  were  critical  to  the  successful  completion  of  this  study.  Mr.  Linker  and  his  staff  at 
the  CBPO  are  acknowledged  for  their  assistance  and  cooperation. 

Mr.  Anthony  S.  Donigian,  Aqua  Terra,  was  the  Principal  Investigator  and  Project  Manager,  responsi¬ 
ble  for  the  overall  technical  direction  of  the  work  and  preparation  of  the  final  report  Mr.  Radha  V. 
Chinnaswamy,  also  of  Aqua  Terra,  reviewed  the  literature  and  assisted  in  developing  nutrient  balances  for 
the  simulated  land-use  categories;  he  also  performed  model  testing  and  calibration  and  assisted  in 
preparing  the  final  report  Messrs.  Brian  BickneU  and  Thomas  Jobes,  Aqua  Terra,  assisted  with 
operational  and  technical  model  application  issues  throughout  the  study. 

At  the  time  of  publication  of  this  report  Director  of  WES  was  Dr.  Robert  W.  Whalin,  and  EL  Director 
was  Dr.  John  Harrison.  Dr.  Richard  E.  Price  was  Chief,  Envirormiental  and  Effects  Division,  EL.  WES 
Commander  was  COL  Robin  R.  Cababa,  EN. 

This  report  should  be  cited  as  follows: 

Donigian,  A.  S.,  Chirmaswaray,  R.  V.,  and  Deliman,  P.  N.  (1998).  “Use  of  nutrient 
balances  in  comprehensive  watershed  water  quality  modeling  of  Chesapeake  Bay,” 

Technical  Report  EL-98-5,  U.S.  Army  Engineer  Waterways  Experiment  Station, 

Vicksburg,  MS. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication. 


VI 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  Units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as  follows: 


^  To  Obtain  Celsius  (C)  temperature  from  Fahrenheit  (F)  readings,  use  the  following 
formula:  C  =  (5/9)(F-32). 


1  Introduction 


Study  Background 

In  1987  the  Chesapeake  Bay  Agreement  was  signed  by  the  EPA  Administrator  and  governors  of  the 
member  states  recommending  a  40%  reduction  in  nutrient  loadings  to  the  Chesapeake  Bay  to  restore 
and  maintain  the  water  quality  and  living  resources  of  the  Bay.  The  Chesapeake  Bay  Watershed 
Model,  based  on  the  U.S.  EPA  Hydrologic  Simulation  Program-Fortran  (HSPF)  (Bicknell  et  al.  1993) 
was  used  to  provide  a  framework  for  quantifying  and  evaluating  the  needed  nutrient  loading 
reductions  to  the  Chesapeake  Bay,  and  to  allow  the  Chesapeake  Bay  Program  Office  (CBPO)  to 
evaluate  the  impacts  of  land  use  changes  and  alternate  nutrient  and  agricultural  management  practices. 
The  Watershed  Model  provides  total  pollutant  loadings  from  all  land  areas  tributary  to  the  Bay  in 
order  to  drive  a  fully  dynamic  three-dimensional,  hydrodynamic/water  quality  model  of  Chesapeake 
Bay. 

The  CBP  Watershed  Model  is  a  unique  state-of-the-art  watershed  modeling  capability  that  includes 
detailed  soil  process  simulation  for  agricultural  areas,  linked  to  an  instream  water  quality  and  nutrient 
model  capable  of  representing  comprehensive  point  and  nonpoint  pollutant  loadings  for  &e  entire 
68,000  square  mile  drainage  area  of  the  Chesapeake  Bay.  In  1992,  major  refinements  and  re¬ 
calibration  of  the  Watershed  Model  was  completed  and  a  report  prepared  describing  the  CBP  Phase  n 
Watershed  Model  application  to  the  Bay  drainage  for  calculating  nutrient  loadings  to  the  Bay 
(Donigian  et  al.  1994).  That  work  included  updating  and  extending  the  model  database,  incorporating 
detailed  agricultural  process  simulation  using  the  AGCHEM  modules  of  HSPF,  developing 
capabilities  for  instream  sediment-nutrient  interactions,  and  re-calibrating  the  improved  model  for  the 
extended  time  period.  However,  for  non-agricultural  areas,  simple  empirically-derived  algorithms 
-were  used  to  estimate  nonpoint  loadings,  both  surface  and  subsurface,  based  on  user-derived  potency 
factors,  washoff  rates,  and  subsurface  concentrations.  This  approach  for  the  non-agricultural  lands, 
including  the  large  areas  of  forested  lands,  did  not  maintain  a  nutrient  balance  and  inhibited  detailed 
assessment  of  the  impacts  and  contributions  of  atmospheric  pollutant  sources  (specifically  nitrogen)  to 
overall  water  quality. 


Scope  and  Objectives 

Recent  HSPF  code  enhancements  (Bicknell  et  al.,  1996a)  have  extended  the  detailed  AGCHEM 
algorithms  to  forested  areas  (Hunsaker  et  al,  1994;  Bicknell  et  al.,  1996b),  have  added  more  direct 
input  of  atmospheric  sources,  and  have  improved  the  AGCHEM  plant  uptake  functions  for  better 
representation  of  agricultural  nutrient  management  praaices  (Donigian  et  al.,  1995).  These 
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enhancements  have  shown  the  added  benefits  of  performing  detailed  nutrient  balance  calculations  for 
the  non-agricultural  areas.  As  a  result  of  these  eiforts,  additional  improvements  and  refinements  were 
identified  and  recommended  to  more  directly  consider  atmospheric  deposition  and  other  nutrient 
sources  by  allowing  nutrient  balance  approaches  for  all  land  uses,  consistent  with  the 
recommendations  of  the  Chesapeake  Bay  Executive  Council  and  the  Nonpoint  Source  Evaluation 
Panel  in  1990  (Chesapeake  Bay  Nonpoint  Source  Evaluation  Panel,  1990).  This  comprehensive 
nutrient  balance  approach  will  help  improve  the  overall  utility  of  the  Watershed  Model  as  a  planning 
tool  for  comprehensive  watershed  planning  and  assessment  of  nutrient  management/reduction 
alternatives. 

The  specific  improvements  recommended  and  tasks  identified  in  this  effort  include  the  following 
developmental  and  application  aspects: 

Model  Development  and  Refinement 

•  Develop  nutrient  balances  and  simulation  procedures  for  the  AGCHEM  module  to  better 
represent  nutrient  cycling,  mass  balance,  and  runoff  contributions  for  non-agricultural 
lands 

•  Test  the  AGCHEM  procedures  and  refinements  on  the  non-agricultural  land  uses,  along 
with  atmospheric  sources,  for  selected  segments  of  the  CBP  Watershed  Model 


Model  Testing  and  Application 

•  Apply  the  refined  AGCHEM  procedures  for  the  non-agricultural  lands  to  a  selected 
subbasin  of  the  Chesapeake  Bay  drainage 

•  Re-calibrate  the  Watershed  Model,  with  the  refined  AGCHEM,  for  the  selected  subbasin 
and  assess  the  load  contributions  from  all  sources  and  the  impact  of  the  refined  procedures 

We  selected  the  Shenandoah  River  subbasin  within  the  Chesapeake  Bay  drainage  to  fully  test  the 
refined  AGCHEM  module  integrated  within  the  Watershed  Model  framework.  Figure  1  shows  the 
CBP  Phase  HI  Watershed  Model  segmentation  for  the  Above  Fall  Line  (AFL)  region,  along  with  the 
location  of  the  Shenandoah  subbasin.  The  most  recent  land  use,  point  source,  meteorologic, 
atmospheric  deposition,  septic  system  load,  flow,  and  water  quality  data  used  in  the  current  Phase  IV 
watershed  model  effort  was  provided  by  the  Chesapeake  Bay  Program  Office  (L.  Linker,  personal 
communication,  1996). 

SUMMARY  CONCLUSIONS  AND  RECOMMENDATIONS 

The  water  quality  results  from  this  effort  are  not  greatly  different  from  those  produced  in  the  earlier 
study,  in  spite  of  all  the  changes  implemented,  including  more  extensive  application  of  AGCHEM, 
addition  of  septic  system  loads,  refinement  of  loading  rates,  and  additional  calibration.  However,  this 
simply  indicates  that  the  prior  calibration  was  a  good  representation  of  the  observed  data,  and  that  the 
changes  implemented  to  better  define  the  load  sources  were  incorporated  while  maintaining  the 
accuracy  of  the  overall  simulation.  The  real  benefits  of  the  current  refinement  phase  of  the  CBP 
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AFL  Model  Segments 


Figure  1 .  Chesapeake  Bay  Watershed  Model  and  Shenandoah  Model  Segments 
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Watershed  Model  are  realized  from  the  extension  of  the  nutrient  balance  approach  to  all  major  land 
use  (except  urban,  in  our  simulations),  and  the  utility  of  this  approach  for  nutrient  management. 

Below  we  discuss  some  of  the  areas  where  the  current  results  differ  from  those  presented  earlier  in 
Donigian  et  al  (1995),  and  identify  some  problems  that  still  remain  for  selected  constituents,  where 
further  ’fine  tuning’  of  the  calibration  is  recommended: 

a.  The  NOj-N  simulation  is  acceptable  but  not  as  seasonally  correct  as  the  earlier  results 
reported  by  Donigian  et  al  (1995).  We  suspect  that  the  differences  primarily  in  late  summer, 
fall,  and  winter  for  selected  years  are  due  mostly  to  the  newly-added  septic  loads  and  possibly 
inaccuracies  in  the  seasonal  loadings  from  the  forest  and  pasture  segments,  i.e.  the  added 
AGCHEM  segments. 

b.  Both  the  NHj-N  and  PO4-P  simulations  have  improved  with  reduced  contributions  from 
cropland  (along  with  reduced  sediment  loads),  reduced  peaks  due  to  nutrient  application 
adjustments,  and  increased  algal  uptake  due  to  higher  benthic  algae  levels  than  in  the  previous 
efforts. 

c.  Simulation  of  both  Organic  N  and  Organic  P  has  improved  in  this  effort,  and  since  Organic  P 
is  the  major  component  of  Total  P  similar  improvement  is  shown  there.  Except  for  the 
November  1985  storm,  the  organics  are  generally  well  represented,  although  many  of  the 
peaks  are  still  somewhat  high.  The  improvements  result  from  reduced  sediment  loads  (and 
associated  organics),  adjustments  to  the  N/P  ratios,  changes  to  manure  applications,  and 
reductions  in  phytoplankton  levels  (which  are  included  in  the  organic  state  variables). 

d.  With  the  detailed  AGCHEM  simulation  of  forests,  the  forest  segment  is  now  a  major 
contributor  of  organic  N,  both  labile  and  refractory  in  both  dissolved  and  particulate  forms, 
along  with  the  animal  acres  segment.  Also,  for  the  forest  segment,  BOD  and  Organic  P  loads 
are  calculated  from  the  organic  N  components;  thus  BOD  loads  from  forest  are  also  a  major 
component.  In  the  Shenandoah,  our  results  indicate  that  forests  contribute  about  25%  of  the 
Organic  N  load  and  30%  of  the  BOD  load,  while  the  animal  acres  segment  contributes  more 
than  30%  of  the  Organic  N  to  the  stream.  Both  of  these  components  need  further 
investigation  as  there  is  little  data  to  confirm  this  level  of  contribution. 


Our  primary  recommendations  resulting  from  this  effort  are  as  follows: 

1.  Finer  segmentation  for  all  stream  reaches  should  be  pursued  as  a  major  component  of  future 
Watershed  Model  enhancements.  Increasing  the  spatial  detail  of  the  model  by  about  a  factor 
of  lOx,  along  with  appropriate  detail  in  the  precipitation  and  land  use  inputs,  will  help  to 
improve  all  the  process  simulations,  with  specific  benefits  for  the  sediment  and  associated 
constituents,  and  benthic  processes. 

2.  A  more  consistent  approach  for  both  BOD  and  organics  loading  needs  to  be  developed  and 
applied  consistently  across  all  land  use  categories.  Currently  the  forest  simulation 
enhancements  provide  loadings  of  both  labile  and  refractory  organic  N  components  (dissolved 
and  particulate),  while  the  other  AGCHEM  segments  (Hi-Till,  Lo-Till,  Hay,  Pasture)  are 
restricted  to  just  the  refractory  particulate  organic  N  and  P  eroded  from  the  land  surface.  The 
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forest  organic  N  capabilities  can,  and  probably  should,  be  applied  to  all  land  segments  to 
implement  this  consistent  loading  representation.  Further  investigation  of  partitioning  and 
transformation  parameters  for  the  organics  will  be  needed,  along  with  consideration  of 
extending  the  forest  N  simulation  approach  to  include  phosphorous  (see  #8,  below). 

3.  The  current  representation  of  septic  system  loads  needs  to  be  re-evaluated.  The  use  of  a 
constant  load  has  helped  to  reduce  the  seasonality  of  the  NO3  simulation  shown  in  both  the 
observed  data  and  the  previous  model  simulations.  A  revised  approach  is  needed  to  allow 
the  septic  loadings  to  be  ’hydrologically-driven’  so  that  the  seasonality  of  the  hydrologic 
regime  and  loadings  is  represented. 

4.  The  algal  simulation,  both  phytoplankton  and  benthic  algae,  need  more  data,  investigation, 
and  evaluation.  In  this  effort,  the  benthic  algal  levels  were  increased  dramatically,  by  factors 
of  10  to  20,  based  on  very  limited  data  from  widely  scattered  sites  outside  the  Chesapeake 
Bay  drainage.  Additional  literature  data  should  be  identified  and  actual  site-specific  data 
within  the  Chesapeake  Bay  watersheds  collected  to  confirm  the  general  magnitude  of  both  the 
benthic  and  phytoplankton  levels  represented  in  the  model.  The  algal  simulation  has  such  a 
critical  impact  on  inorganic  nutrient  levels  that  major  improvements  in  their  modeling  will 
depend  on  establishing  realistic  levels  for  the  algal  populations. 

5.  The  urban  land  use  should  be  divided  into  separate  urban  categories  —  residential, 
commercial,  industrial  —  each  with  pervious  and  impervious  fractions,  as  a  prelude  for 
nutrient  mass  balance  and  AGCHEM-type  model  application.  Although  CBPO  has  pursued  an 
AGCHEM  approach  for  the  aggregated  urban  pervious  segment,  we  feel  a  better  definition  of 
specific  urban  activities  is  needed  to  develop  reasonable  nutrient  balances. 

6.  The  CBPO  should  explore  the  option  of  eliminating  the  current  ’composite  crop’ 
representation  in  the  model,  developed  as  part  of  the  Phase  11  enhancements  in  1991,  and 
move  to  simulating  each  major  crop  individually  in  each  cropland  category.  Recent  computer 
hardware  developments  have  eliminated  many  of  the  run  time  restrictions  that  required  this 
simplification,  and  such  an  approach  would  allow  more  accurate  representation  of  agricultural 
practices  and  the  resulting  nutrient  balances. 

7.  In  conjunction  with  representing  each  major  crop,  the  nutrient  application  rates,  timing, 
procedures,  and  composition  distribution  (both  fertilizers  and  manure)  should  be  closely 
reviewed  and  revised  as  needed.  Experience  with  both  the  Watershed  Model  and  the  detailed 
Patuxent  Model  has  confirmed  the  critical  importance  of  the  assumptions  underlying  the 
nutrient  applications  in  the  model. 

8.  The  forest  N  simulation  approach  should  be  extended  to  include  the  P  cycle,  so  that  both  N 
and  P  mass  balances  can  be  implemented  for  all  land  segments.  Just  as  the  previous 
AGCHEM  module  provided  a  valid  framework  for  the  forest  N  enhancements,  the  P  cycle 
processes  currently  in  AGCHEM  can  be  readily  adapted  for  forested  conditions.  Field  site 
testing  on  small  forested  watersheds  would  be  need^  to  fully  evaluate  the  code  enhancements. 


Chapter  1  Introduction 


Format  of  this  Report 


Following  this  introduction,  Chapter  2  describes  the  development  of  nutrient  balances  for  all  land  uses 
with  the  focus  on  the  non-agricultural  cropland  segments  while  Chapter  3  presents  the  results  of 
model  testing  and  re-calibration  to  the  Shenandoah  subbasin.  The  Appendix  includes  complete 
simulation  results  for  the  Shenandoah  subbasin.. 
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Chapter  1  Introduction 


2  Expected  Nutrient  Balances  by  Land 
Use  Categories 


Overview  and  Summary  Tables 

As  part  of  this  study,  ’expected’  nutrient  balances  were  developed  for  forest,  pasture,  and  urban  land 
use  categories  to  help  guide  the  model  application  and  calibration  effort  for  those  land  segments  for 
which  the  new  HSPF  Version  11.0  AGCHEM  module  was  applied  to  simulate  the  detailed  nitrogen 
dynamics.  Consequently,  the  AGCHEM  sections  replaced  the  PQUAL  sections  in  the  Watershed 
Model  for  forest  for  only  the  N  species,  and  for  pasture  for  both  N  and  P  species.  Although  the 
preliminary  nutrient  balances  for  the  urban  land  use  were  developed,  the  nutrient  loadings  from  the 
urban  segment  in  the  Model  were  still  simulated  using  PQUAL,  as  in  Phase  HI;  the  basis  for  this 
decision  is  discussed  in  Chapter  3. 

Tables  1  and  2  (at  the  end  of  this  chapter)  show  the  typical  or  expected  nitrogen  and  phosphorus 
balances,  respectively,  for  different  major  crops  and  land  cover/use  categories.  The  information  in 
the  tables  is  presented  in  a  ’production’  sense  by  estimating  the  annual  INPUTS  and  OUTPUTS  for 
the  soil-plant  system.  The  INPUTS  represent  external  additions  to  the  system,  such  as  nutrient 
applications  (i.e.  fertilizer  and  manure),  in  addition  to  net  mineralization  from  the  soil  that  supplies 
plant-available  inorganic  nutrients.  The  OUTPUTS  represent  various  loss  mechanisms,  plus  plant 
uptake  (e.g.  through  harvest  or  plant  retention)  that  extracts  the  nutrients  from  the  soil  impacting  the 
potential  for  nutrient  export  and  losses.  Thus,  although  mineralization  and  plant  uptake  are  not  truly 
external  to  the  soil-plant  system,  they  are  most  often  key  components  in  establishing  representative 
nutrient  balances  for  most  land  use/cover  categories. 

A  review  of  the  literature  was  performed  to  develop  the  nutrient  balances  for  each  category.  The 
nitrogen  balance  for  forest  was  derived  primarily  from  information  in  the  report  by  Oak  Ridge 
National  Laboratory  (ORNL;  Hunsaker  et  al.,  1994)  that  was  the  basis  for  the  algorithm 
enhancements  to  AGCHEM  for  the  forest  N  cycling  (Bicknell  et  al.,  1996).  Since  HSPF  does  not 
currently  have  the  capability  to  simulate  the  forest  P  cycling  in  detail,  the  PQUAL  module  was  used 
for  P  loading  simulations  and  hence  we  did  not  develop  the  P  balance  for  forests.  Also,  the  nutrient 
balances  for  agricultural  croplands  developed  during  Phase  II  (Donigian  et  al.,  1994)  were  used  in 
this  study  and  are  included  in  Tables  1  and  2.  In  this  chapter  (below),  we  briefly  discuss  the 
development  of  the  ’expected’  nutrient  balances  for  each  land  use/cover  as  background  for  the  values 
shown  in  the  tables. 


Chapter  2  Expected  Nutrient  Balances 
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Forest  Balances 


In  the  previous  CBP  modeling  efforts,  the  PQUAL  module  in  HSPF  was  used  to  compute  the  N  and 
P  loadings  from  the  forest  segments  based  on  user-defined  potency  factors  and  subsurface 
concentrations.  Since  the  Chesapeake  Bay  watershed  has  approximately  60%  forested  land  and  there 
was  concern  regarding  the  level  of  N  loadings  from  these  areas,  there  was  a  need  to  simulate  the  N 
cycling  in  forests  with  a  more  detailed  mass  balance  approach  that  would  also  allow  more  direct 
consideration  of  atmospheric  deposition  inputs  and  impacts. 

The  ORNL  report  (Hunsaker  et  al.,  1994)  provided  the  design  details  for  the  enhancement  of  the 
AGCHEM/NTTR  module  in  HSPF  based  on  an  extensive  literature  review  of  forest  nitrogen  pools 
and  fluxes,  review  of  monitored  N  data  from  the  CBP  region,  and  available  N  models.  The  literature 
review  focused  on  the  data  collected  from  two  studies;  the  International  Biological  Program  which 
provided  data  on  116  forest  research  sites  around  the  world,  and  the  Integrated  Forest  Study  that 
provided  17  forest  research  sites  (16  in  North  America).  The  information  and  data  collected  at  these 
sites  and  subsequently  presented  in  the  ORNL  report  formed  a  reasonable  basis  for  developing  the 
nitrogen  balance  shown  in  Table  1,  besides  serving  as  a  tool  for  the  development  of  forest  N  module. 
Based  on  the  review  of  available  literature  and  information  provided  in  ORNL  report,  the  following 
general  summary  is  presented. 

•  The  primary  sources  of  input  to  forests  are  atmospheric  N  deposition,  with  nitrogen  fixation 
for  some  species,  along  with  plant  available  N  (inorganic  N)  from  mineralization.  Forest 
fertilization  can  be  important  in  silvicultural  activities  and  should  be  included  if  appropriate 
for  a  specific  site  assessment;  we  have  not  included  forest  fertilization  in  Table  1. 

•  The  major  pathways  by  which  N  export  losses  occur  include  leaching  from  soil  and 
denitrification.  Also,  the  surface  runoff  losses  are  generally  small  except  when  the  forest 
system  has  reached  higher  levels  of  N  saturation  (Stoddard,  1994).  Although  plant  uptake  is 
not  a  true  ’loss’  from  the  system,  it  is  a  key  component  in  the  overall  balance. 

•  The  other  processes  that  play  important  roles  in  forest  N  cycling  are  retention  of  ammonium- 
N  by  soils,  immobilization  of  available  nitrogen  by  microorganisms,  and  return  of  plant  N  to 
the  soil  both  belowground  and  through  the  forest  litter  layer.  See  Hunsaker  et  al  (1994)  and 
Bicknell  et  al  (1996c)  for  additional  details. 

•  Since  mineralization  and  plant  uptake  are  such  dominant  components  of  a  forest  N  balance,  an 
accurate  accounting  of  their  fluxes  is  critical  to  modeling  the  N  cycling  and  export  to 
waterbodies. 


Pasture  Balances 

A  significant  portion  of  the  Chesapeake  Bay  region  consists  of  pasture,  or  grassland  ecosystems.  The 
grassland/pasture  ecosystem  presents  a  wide  diversity  of  environments,  productivity  and  degree  of 
management.  The  limited  literature  data  available  for  the  input-output  dynamics  of  nutrients  in  both 
grasslands  and  pastures  are  often  based  on  studies  conducted  across  the  continental  U.S.  and  are  not 
specific  to  CBP  region.  However,  based  on  the  review  of  available  literature  data  on  nitrogen  cycling 


8 


Chapter  2  Expected  Nutrient  Balances 


in  these  ecosystems,  the  following  general  conclusions  are  presented.  Since  there  are  very  little 
literature  data  available  on  phosphorus  cycling  in  the  grassland  systems,  ,many  of  the  corresponding 
components  were  assumed  to  be  analogous  to  N  cycle  components. 

•  The  principal  sources  of  N  inputs  are  atmospheric  deposition,  manure  and  N  fertilizer 
applications,  along  with  the  mineralization  contribution  to  plant-available  N  as  noted  above  for 
forests.  Symbiotic  and  nonsymbiotic  fixation  of  N2  is  considered  small  or  insignificant  under 
typical  grassland  conditions  (Woodmansee,  1978). 

•  The  N  uptake  by  plants  (aboveground,  belowground  and  understory)  are  in  the  range  of  65- 
80%  of  the  total  N  input  (Legg  and  Meisinger,  1982;  Muchovej  and  Rechcigl,  1994).  Hence, 
an  average  73%  plant  N  uptake  was  assumed  in  Table  1. 

•  In  pasture,  loss  due  to  volatilization  is  significant  (Muchovej  and  Rechcigl,  1994).  About  10- 
25%  is  lost  due  to  volatilization  of  NH3  from  animal  urine  and  feces  while  about  5%  is  lost 
due  to  surface  runoff  (Meisinger  and  Randall,  1991;  Legg  and  Meisinger,  1982). 

Consequently,  an  average  17%  volatilization  loss  and  5%  surface  runoff  loss  were  assumed  in 
Table  1;  clearly,  these  values  will  vary  by  site-specific  conditions. 

•  Denitrification  is  not  a  significant  factor  in  N  balance  for  most  pasture  and  grasslands 
(Woodmansee,  1978).  This  is  probably  due  to  the  fact  that  native  and  extensively  managed 
grasslands  are  N  deficient  (Muchovej  and  Rechcigl,  1994).  However,  denitrification  is  very 
site-specific  and  was  assumed  to  be  similar  to  hay  land  for  our  balance. 

•  Some  studies  indicate  that  leaching  is  an  important  factor  while  other  studies  indicate  leaching 
is  either  small  or  insignificant  in  these  systems  (Woodmansee,  1978;  Legg  and  Meisinger, 
1982;  Keeney,  1982).  Since  very  little  research  has  been  done  on  nitrate  leaching  from 
grasslands  in  the  U.S.  (Muchovej  and  Rechcigl,  1994)  and  those  studies  were  conducted  at 
different  soils,  climatic  conditions  and  vegetation  type,  it  is  difficult  to  generalize  the  leaching 
potential  of  pasture  land.  However,  subsurface  losses  are  one  of  the  major  pathways  by 
which  N  is  transported  from  grasslands  (Muchovej  and  Rechcigl,  1994).  Hence,  the  leaching 
and  subsurface  losses  were  assumed  to  be  the  same  as  that  for  hay  lands  since  hay  and  pasture 
lands  have  similar  conditions. 

•  Mineralization  and  inunobilization  are  important  factors  in  the  grassland  N  balance.  However, 
data  are  not  available  on  the  mineralization  and  immobilization  rates  for  grassland  systems. 

As  mentioned  above,  since  pasture  and  hay  lands  have  similar  conditions,  mineralization  and 
immobilization  for  pasture  were  assumed  to  be  in  the  same  range  as  that  for  hay  lands.  For 
conditions  where  manure  inputs  are  significant,  mineralization  rates  in  the  surface  soil  layers 
should  reflect  the  higher  values  common  to  the  readily  mineralizable  portion  of  manure 
organics. 


Urban  Balances 

There  are  very  few  studies  and  literature  data  on  the  nutrient  dynamics  in  an  urban  environment. 
Although  the  few  available  studies  were  performed  in  parts  of  the  northeastern  U.S.,  they  provide 
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current  state  of  knowledge  on  the  fate  of  the  fertilizers  applied  to  home  lawns.  Based  on  the  review 
of  available  literature  data,  the  following  general  conclusions  are  presented.  Since  there  are  no 
literature  data  available  on  phosphorus  cycling  in  urban  land  use,  it  was  assumed  to  be  analogous  to 
N  cycling  in  our  tables. 

•  The  principal  source  of  N  input  is  fertilizer  N  applied  to  lawns.  Intensively  managed 
turfgrass  receives  between  100-200  lb  N/ac/yr  (Muchovej  and  RechcigI,  1994;  Petrovic, 

1990;  Morton  et  al.  1988;  Gold  et  al.,  1990). 

•  In  general,  the  N  uptake  by  turfgrass  is  in  the  range  of  5-74%  of  the  total  N  input  depending 
on  the  N  release  rate,  application  rate,  and  species  of  grass  (Petrovic,  1990).  On  an  average, 
about  35-60%  of  the  applied  N  is  found  in  clippings  while  about  14-21%  of  the  fertilizer  N  is 
found  in  a  thatch  layer  (Petrovic,  1990;  Muchovej  and  RechcigI,  1994).  Therefore,  the  total 
N  taken  up  by  turfgrass  is  in  the  range  of  50-80%.  In  Table  1,  we  assumed  that  an  average 
65%  of  the  total  N  input  was  taken  up  by  plants. 

•  Lx)ss  of  applied  fertilizer  to  the  atmosphere  as  either  ammonia  due  to  volatilization  or  as  one 
of  several  nitrous  oxide  compounds  (e.g.  denitrification)  is  a  significant  factor  in  the  N 
balance  of  turfgrass  (Petrovic,  1990).  Factors  such  as  presence  or  absence  of  thatch, 
irrigation  and  humidity  can  affect  the  rate  of  volatilization.  The  soil  moisture  and  temperature 
have  a  significant  influence  on  denitrification  rates.  About  10-36%  of  the  applied  fertilizer  N 
is  lost  due  to  volatilization  and/or  denitrification  (Petrovic,  1990).  In  Table  1,  we  assumed 
that  an  average  23%  was  lost  due  to  volatilization  and  denitrification. 

•  The  fertilizer  management  practices,  soil  texture,  and  irrigation  appear  to  have  influence  on 
the  leaching  losses  in  turfgrass.  Even  though  the  leaching  losses  from  turfgrasses  are 
variable,  in  general  these  losses  are  less  than  10%  of  the  applied  N  (Muchovej  and  RechcigI, 
1994;  Morton  et  al.  1988).  Hence,  we  assumed  that  ten  percent  of  the  total  N  input  is  lost 
due  to  leaching. 

•  Surface  runoff  losses  are  minimal  (less  than  7  %  of  the  total  waterborne  loss)  in  turfgrass  with 
permeable  soils  (Morton  et  al.  1988;  Petrovic,  1990).  For  our  purposes,  we  assum^  that 
about  5%  of  the  total  N  input  is  lost  in  surface  runoff. 

•  Mineralization  rates  are  not  known  or  little  research  has  been  done  for  turfgrass  (Petrovic, 
1990).  However,  we  assumed  that  N  and  P  mineralization  rates  for  urban  lawns  are  in  the 
same  range  as  that  for  hay  and  pasture  land.  Studies  comparing  rural  and  urban  forest  N 
cycling  indicate  lower  mineralization  rates,  less  cycling,  and  lower  labile  N  and  nitrification 
in  the  urban  environments  possibly  due  to  pollution  effects  on  vegetation  and  litter  quality 
(Goldman  et  al.,  1995;  White  and  McDonnell,  1988). 

A  significant  limitation  in  applying  the  urban  nutrient  balance  within  the  framework  of  the  Watershed 
Model  is  the  lumped,  or  aggregate,  nature  of  the  urban  land  use  segment  that  includes  residential, 
commercial,  industrial,  parks,  etc.  A  more  refined  definition  of  the  urban  category,  with  the  specific 
activities  and  land  cover  is  needed  to  attempt  a  nutrient  balance  modeling  approach  (see  Chapter  3). 
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Agricultural  Cropland  Balances 


The  nitrogen  and  phosphorus  balances  developed  for  croplands  during  the  Phase  II  modeling  effort 
(Donigian  et  al.,  1994)  were  also  used  in  this  study  without  modification.  The  values  included  in 
Tables  1  and  2  are  typical  nitrogen  and  phosphorus  balances  expected  for  Com,  Soybean,  Grains  and 
Hay  when  the  nutrients  are  applied  at  agronomic  rates  meeting  crop  requirements.  Moreover,  since 
plant  uptake  amounts  are  a  function  of  crop  yields,  die  uptake  values  are  based  on  average  yields 
expected  for  the  Chesapeake  Bay  region.  Our  experience  with  using  these  values  for  agricultural 
croplands  indicates  that  wide  variations  can  be  expected  depending  on  site-specific  conditions.  For 
example,  yields  can  vary  with  species,  double-cropping  and  winter  cover  crops  will  increase  uptake 
levels,  and  the  composition  of  nutrient  inputs  (e.g.  organic  versus  inorganic),  especially  for  manure 
and  sludge  applications,  will  have  a  major  impact  on  the  N  and  P  balance  components.  Thus,  users 
should  attempt  to  develop  site-specific  balances  to  the  extent  possible  and/or  refine  the  values  in  the 
tables  for  specific  local  conditions. 


Closure 

The  nutrient  balances  in  Tables  1  and  2  for  forest,  pasture  and  urban  land  uses  were  based  on 
literature  data  collected  in  different  parts  of  the  U.S.,  but  with  a  specific  focus  on  the  Chesapeake  Bay 
region,  while  the  values  for  the  croplands  are  more  directly  pertinent  to  the  Bay  region..  Even 
though  the  data  may  have  been  derived  from  somewhat  diverse  sites  and  conditions,  and  the 
techniques  used  in  collecting  the  data  may  have  been  different,  the  available  literature  and  the  data 
formed  a  reasonably  good  basis  for  determining  the  general  magnitude  of  many  of  the  nutrient 
balance  components.  Model  users  should  use  these  balances  as  general  guides  in  evaluating 
simulation  results  with  the  knowledge  that  local  conditions  can  have  a  major  impact  on  the  values 
shown  in  the  tables.  As  noted  above,  these  balances  should  be  updated  and/or  modified  whenever 
more  site-specific  data  are  available. 


Chapter  2  Expected  Nutrient  Balances 


11 


Table  1.  Typical  Nitrogen  Balance  For  Major  Crops  and  Land  Use/Land  Cover  Categories  (Ib/ac/yr) 

Corn*  Soybeans*  Grains*  Hay*  Forest’ 
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Table  2.  Typical  Phosphorus  Balances  For  Major  Crops  and  Land  Use/Land  Cover  Categories  (Ib/ac/yr) 
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3  Application  of  Nutrient  Balances  and 
Model  Refinements  to  the  Shenandoah 
River  Watershed 


Overview  of  Model  Refinements 

The  model  testing  and  refinement  focussed  on  the  evaluation  of  the  new  forest  simulation  procedures 
and  AGCHEM  application  to  the  non-cropland  segments,  but  also  included  other  enhancements  to  the 
modeling  that  were  included  as  part  of  the  refinements  under  Phase  IV  of  the  CBP  Watershed  Model 
effort.  The  Phase  III  model  refinement  effort,  discussed  and  reported  by  Donigian  et  al  (1995)  for 
three  selected  test  basins,  included  the  following  changes  from  the  earlier  Phase  II  work  (Donigian  et 
al,  1994): 

a.  Land  use  was  updated  by  CBPO  to  1990  conditions  (Neumiller  at  al.  1994). 

b.  Selected  model  segments  were  re-defined  for  finer  spatial  detail. 

c.  Precipitation  and  meteorologic  data  was  reviewed  and  extended  through  1991. 

d.  Point  sources,  diversion,  and  atmospheric  deposition  files  were  extended  through  1991. 

e.  The  yield-based  plant  uptake  function  was  used  in  AGCHEM. 

f.  Enhanced  SPECIAL  ACTIONS  capabilities  were  used  to  reduce  the  length  of  the  model  input 
and  improve  representation  of  chemical  application  practices. 

g.  Atmospheric  deposition  was  included  as  timeseries  directly  input  to  chemical  storages  for  the 
AGCHEM  segments,  taking  advantage  of  additional  HSPF  Version  No.  11  capabilities. 

h.  Benthic  oxygen  demand  and  benthic  algae  processes  were  activated  within  HSPF  in  order  to 
improve  the  low  DO  and  inorganic  nutrient  simulation. 

As  part  of  this  study,  in  coordination  with  the  Phase  IV  efforts  of  the  CBPO  staff,  the  following 
additional  changes  and  model  refinements  were  implemented: 
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Forest  nitrogen  simulation  capabilities  (implemented  within  AGCHEM)  in  HSPF  Version  No. 
11  (Bicknell  et  al.,  1996a)  were  used  for  die  forest  segment..  The  forest  phosphorous 
simulation  is  the  same  as  in  Phase  III,  i.e.  the  PQUAL  module  of  HSPF  is  us^  for 
phosphorous  loadings. 

AGCHEM  is  used  for  both  N  and  P  simulations  for  pasture,  in  addition  to  the  cropland  and 
hay  segments. 

Atmospheric  deposition  was  added  for  all  segments,  including  urban  impervious. 

Septic  system  loads  of  N  were  estimated  and  input  as  a  point  load  for  each  stream  reach 
segment,  as  recommended  by  CBPO  (L.  Linker,  personal  communication,  1996). 

Minor  adjustments  were  made  to  the  SPECIAL  ACTIONS  operations  to  improve  the 
representation  of  the  nutrient  applications  to  croplands,  hay,  and  pasture. 

Manure  nutrient  applications  in  Phase  II/III  were  reviewed  by  CBPO  based  on  a  manure  mass 
balance  analysis;  application  rates  were  revised  and  values  for  pasture  were  developed  (M. 
Palace,  personal  communication,  1996). 

Water  quality  calibration  included  a  greater  focus  on  the  impact  of  benthic  algae  processes  on 
inorganic  nutrient  concentrations,  in  addition  to  review  and  assessment  of  all  water  quality 
constituents. 

The  Watershed  Model  re-calibration  and  refinement  effort  for  the  Shenandoah  segments  involved: 

a.  Review,  calibration,  and  evaluation  of  the  simulations  for  the  new  AGCHEM  applications  to 
the  forest  and  pasture  segments,  and  a  check  of  the  cropland  simulations,  based  on  the 
expected  nutrient  balances  developed  in  this  study 

b.  Review  and  comparison  of  nonpoint  source,  point  source,  atmospheric  deposition,  and  septic 
loadings  from  all  land  use  categories 

c.  Calibration  of  selected  instream  water  quality  parameters  based  on  comparison  of  simulated 
and  observed  concentrations 

These  tasks  and  simulation  results  are  discussed  below  under  the  ’Nonpoint  Source  and  Loading 
Assessment’  and  ’Water  Quality  Calibration  Results’  sections,  following  a  brief  discussion  of  the 
’Hydrology  Simulation’  results.  This  chapter  concludes  with  our  ’Conclusions  and 
Recommendations’  derived  from  the  results  of  these  efforts  and  our  interactions  with  the  CBPO  staff 
in  their  application  of  the  Phase  IV  refinements  to  the  entire  Chesapeake  Bay  Watershed. 


Hydrology  Simulation 

The  hydrology  simulation  for  this  study  on  the  Shenandoah  Basin  is  essentially  the  same  as  was 
reported  by  Donigian  et  al.,  (1995)  for  the  prior  study,  with  the  following  exceptions: 
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a.  The  CBPO  staff  reviewed  the  precipitation  data  for  the  entire  Chesapeake  Bay  Watershed 
Model  and  updated  selected  datasets  based  on  an  analysis  of  spatial  variation  and  changes  with 
elevation  for  major  storm  events  (Neumiller  et  al.,  1994).  Except  for  a  few  storm  events, 
these  adjustments  did  not  have  a  major  impact  on  the  hydrology  simulation. 

b.  As  part  of  the  forest  nitrogen  enhancements  to  HSPF  (Bicknell  et  al.,  1996c),  model  testing 
was  performed  on  completely  forested  sites  resulting  in  a  refinement  of  the  hydrology 
parameters  for  the  forest  segments.  These  refinements  only  affected  the  distribution  between 
surface  and  subsurface  flow,  with  no  significant  impact  on  total  annual  flow  volumes 
(Donigian  and  Chinnaswamy,  1996).  The  new  parameters  produced  less  surface  runoff  and 
more  interflow  than  previously  simulated,  providing  a  better  representation  of  the  hydrologic 
regime  for  forested  watersheds. 

c.  In  this  study,  we  reviewed  the  hydrology  simulation  at  two  sites  within  the  Shenandoah  Basin; 
the  Basin  outlet  at  Millville,  WV,  and  the  South  Fork  Shenandoah  River  at  Front  Royal,  VA. 
In  addition  to  the  Shenandoah  River  at  Millville,  which  has  been  the  calibration  station  for  the 
previous  modeling  efforts,  we  output  model  results  for  the  Reach  190  which  corresponds  to 
the  Front  Royal  station.  No  calibration  was  needed,  as  shown  by  the  results  below,  as  the 
model  parameters  used  for  the  entire  Basin  were  appropriate  for  the  upper  model  segment. 

Table  3  shows  the  comparison  of  simulated  and  observed  flow  volumes  for  each  year  and  for  the 
entire  simulation  period  for  both  sites  within  the  Shenandoah  Basin.  Figures  2  and  3  show  the  daily 
flow  simulation  results,  and  Figures  4  and  5  show  the  flow  duration/frequency  comparisons  for  both 
sites,  respectively.  The  Appendix  also  includes  these  results  along  with  the  daily  flow  comparison  at 
an  expanded  scale  for  a  clearer  definition  of  the  two  curves. 

The  hydrology  simulation  results  are  clearly  a  strength  of  the  CBP  Watershed  Model.  The  annual 
volumes,  flow  frequency,  and  daily  timeseries  comparisons  all  show  very  good  agreement  between 
simulated  and  observed  values.  As  in  any  modeling  exercise,  some  differences  are  evident  in  each 
comparison,  but  the  overall  agreement  is  very  good  and  provides  a  sound  basis  for  the  water  quality 
simulation  which  is  the  focus  of  this  effort. 

Nonpoint  Source  and  Loading  Assessment 

The  loading  assessment  involved  reviewing  the  contaminant  loadings  for  each  constituent  for  each 
source,  including  nonpoint  sources,  point  sources,  septic  loads,  and  atmospheric  deposition.  Tables  4 
and  5  are  key  tools  in  this  assessment;  Table  4  shows  the  unit  area  (per  acre)  loads  for  each  land  use 
simulated  in  the  Watershed  Model  for  Segments  190  and  200  within  the  Shenandoah  Basin,  while 
Table  5  shows  the  ’percent  of  total  load’  derived  from  each  source.  These  tables  show  only  the  mean 
values  for  the  eight-year  simulation  period,  with  the  complete  results  provided  in  the  Appendix;  the 
Appendix  also  includes  complete  summaries  of  the  AGCHEM  simulations  for  the  forest,  cropland 
(Hi-Till,  Lo-Till,  Hay),  and  pasture  segments. 


16 


Chapter  3  Application  to  the  Shenandoah  River  Watershed 


Table  3.  Shenandoah  River  Watershed  Hydrologic  Calibration:  Comparison  of  Annual 
Total  Observed  vs  Simulated  Flow 


SF  SHENANDOAH  RIVER  AT  FRONT  ROYAL,  VA  (SEGMENT  190) 


YEAR 

OBSERVED 
FLOW  (in) 

SIMULATED 
FLOW  (in) 

1984 

18.60 

20.05 

1985 

16.20 

17.10 

1986 

7.34 

7.24 

1987 

16.10 

14.73 

1988 

7.71 

7.06 

1989 

15.08 

16.35 

1990 

13.58 

15.41 

1991 

11.67 

11.63 

MEAN 

13.29 

13.70 

SHENANDOAH  RIVER  AT  MILLVILLE,  WV  (SEGMENT 

YEAR 

OBSERVED 

SIMULATED 

FLOW  (in) 

FLOW  (in) 

1984 

17.82 

18.27 

1985 

12.99 

14.35 

1986 

6.92 

6.88 

1987 

13.95 

12.09 

1988 

7.72 

6.72 

1989 

12.02 

13.15 

1990 

11.58 

13.78 

1991 

10.50 

9.99 

MEAN 

11.69 

11.90 
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FIGURE  +.  FREQUENCY  ANALYSIS  OF  FLOW  AT  REACH  190 
SF  SHENANDOAH  RIVER  AT  FRONT  ROYAL.  VA 


Table  4  allows  the  model  user  to  assess  validity  and  reasonability  of  the  unit  area  loading  values  by 
comparing  them  to  ’expected’  or  generally  accepted  values,  or  ranges  of  values,  from  the  literature. 
Table  5  allows  the  user  to  determine  the  relative  distribution  and  contributions  of  each  load  source  to 
each  stream  reach  in  the  Watershed  Model.  Judicious  and  careful  analysis  of  these  two  types  of 
tables  for  each  model  segment  allows  the  user  to  evaluate  both  the  validity  of  the  loading  sources  and 
identify  those  sources  that  may  need  to  be  re-evaluated  in  order  to  improve  or  refine  the  comparison 
between  simulated  and  observed  instream  concentrations  which  are  a  direct  result  of  these  loads. 

In  this  study,  the  nonpoint  source  and  loading  assessment  involved  application  of  a  nutrient  mass 
balance  modeling  approach  to  forest  and  pasture  land  segments,  a  review  and  relative  comparison  of 
the  nonpoint  loading  rates  for  all  the  land  uses,  and  an  assessment  of  the  relative  load  contributions 
from  all  sources  as  derived  from  the  information  in  Table  5.  These  efforts  involved  the  following 
components: 

a.  As  discussed  in  Chapter  2,  the  literature  was  reviewed  to  develop  ’expected’  nutrient  balances 
for  all  land  uses,  with  a  primary  focus  on  forest,  pasture,  and  urban,  as  these  categories  had 
not  been  modeled  with  an  AGCHEM,  nutrient  balance  approach  in  prior  phases  of  the  CBP 
Watershed  Model.  The  N  and  P  balances  shown  in  Tables  1  and  2  of  Chapter  2  provided  a 
general  target  for  the  nonpoint  simulation  of  each  land  use,  recognizing  that  significant 
deviations  are  possible  due  to  local  and  regional  conditions. 

b.  The  forest  N  modeling  with  the  enhanced  AGCHEM  module  was  derived  from  test 
applications  of  the  model  enhancements  to  selected  small  forested  watersheds  in  and  near  the 
Shenandoah  Basin,  and  a  regional-scale  application  as  part  of  the  code  development  and 
testing  effort  (Bicknell  et  al.,  1996c).  One  of  the  key  conclusions  of  the  Phase  II  loading 
rate  assessment  was  the  identification  of  the  relatively  high  loading  rates  of  NOj  from  the 
forest  segment.  In  comparison  with  selected  datasets  for  forested  watershed  sites,  Hunsaker, 
Garten,  and  Mulholland  (1994)  noted  that  the  NO3  loading  rates  used  in  the  Watershed  Model 
were  up  to  ten  times  greater  in  the  Susquehanna  than  indicated  by  measured  data,  and  up  to 
five  times  greater  in  other  model  segments.  The  model  values  were  typically  in  the  range  of 
1.0  to  6.0  Ib/ac/year,  whereas  the  measured  data  were  in  the  range  of  0.5  to  2.0  Ib/ac/yr. 
Therefore  in  this  effort  we  adjusted  the  forest  AGCHEM  parameters,  within  the  overall 
expected  N  balance,  to  reduce  the  model  generated  forest  NO3  loads  to  the  general  range  of 
1.0  to  2.0  Ib/ac/yr. 

c.  The  pasture  segment  was  converted  to  nutrient  balance  modeling,  using  AGCHEM  for  both  N 
and  P,  from  the  empirical  PQUAL  approach  used  in  earlier  phases  of  the  Watershed  Model. 
Pasture  segment  parameters  were  initially  derived  from  the  Hay  segment,  with  application 
rates  and  procedures  for  N  provided  by  the  CBPO  (M.  Palace,  CBPO,  personal 
communication,  1996).  We  developed  corresponding  P  application  rates  based  on  the  N/P 
ratio  used  in  the  Phase  II/III  manure  nutrient  application  rates.  The  yield-based  plant  uptake 
option  was  used  and  calibration  was  performed,  primarily  on  the  N  and  P  uptake  targets,  and 
within  the  general  guidance  of  the  expected  nutrient  balances  for  pasture. 

d.  The  data  and  information  on  nutrient  balances  for  urban  land  uses  is  extremely  limited  and 
anecdotal.  Furthermore,  since  our  urban  category  is  an  aggregation  of  all  relevant  urban  land 
uses  within  a  model  segment,  developing  an  appropriate  nutrient  balance  was  problematic 
without  a  clear  definition  of  the  specific  activities  taking  place  on  the  urban  segment. 
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Table  4.  Unit  Area  Nonpoint  Source  Loading  Rates  for  Each  Land  Use  for  the  Shenandoah 
Basin  (Ib/ac,  except  SED  which  is  tons/ac) 


CONSTITUENT 


< - 

FOR 


- Pervious - . - 

HTC  LTC  PAS  URB 


— >< — Impervious — >  Total 
HAY  ANML  RES  Load 


Segment  190 


Segment  200 


NH3 

0.078 

2.726 

2.323 

0.544 

0.308 

0.495 

221.504 

2.145 

1.669 

NOS 

0.880 

17.240 

14.892 

4.340 

5.788 

3.997 

55.376 

6.291 

4.394 

ORGN 

1.020 

3.324 

2.842 

1.303 

1.972 

0.683 

1661.278 

3.293 

2.113 

TN 

1.978 

23.289 

20.055 

6.188 

8.068 

5.175 

1605.902 

11.729 

8.041 

P04 

0.022 

1.735 

1.475 

0.230 

0.255 

0.395 

55.376 

0.429 

0.538 

ORGP 

0.025 

0.927 

0.762 

0.245 

0.282 

0.186 

83.064 

0.470 

0.263 

TP 

0.047 

2.662 

2.237 

0.476 

0.537 

0.581 

276.880 

0.899 

0.857 

BOD 

9.995 

86.490 

51.653 

6.796 

15.949 

14.243 

3876.315 

26.625 

17.403 

SED 

0.057 

0.979 

0.676 

0.182 

0.152 

0.253 

0.000 

0.000 

0.165 

NH3 

0.057 

2.007 

1.632 

0.353 

0.212 

0.317 

185.347 

1.949 

1.244 

N03 

1.042 

13.186 

11.305 

3.281 

4.598 

3.255 

46.337 

5.939 

3.514 

ORGN 

0.742 

2.593 

1.850 

1.046 

1.276 

0.249 

1390.104 

3.236 

1.403 

TN 

1.841 

17.785 

14.786 

4.680 

6.085 

3.821 

1343.767 

11.124 

6.071 

P04 

0.005 

1.345 

1.153 

0.135 

0.162 

0.269 

46.337 

0.421 

0.288 

ORGP 

0.018 

0.726 

0.496 

0.195 

0.182 

0.067 

69.505 

0.462 

0.162 

TP 

0.023 

2.071 

1.648 

0.330 

0.344 

0.336 

231 .684 

0.883 

0.488 

BOD 

7.265 

61.286 

29.006 

5.558 

10.320 

8.100 

3243.574 

26.168 

11.795 

SED 

0.039 

0.770 

0.441 

0.140 

0.143 

0.095 

0.000 

0.000 

0.108 

Table  5.  Percent  OF  Total  Load  Contributed  From  Each  Source  in  Shenandoah  Basin 


con: 

STITI 

JENT  <• 

---Perv 

lOUS---- 

•Impervious-> 

Atmos 

Point 

Septic 

Total 

FOR 

HTC 

LTC 

PAS 

URB 

HAY 

ANML 

RES 

Dep 

Source 

Load 

Load 

Segment 

190 

NH3 

2.33 

3.60 

7.13 

7.12 

1.32 

3.55 

5.39 

2.69 

0.39 

66.47 

0.00 

100.00 

N03 

9.93 

8.66 

17.37 

21.58 

9.41 

10.88 

0.51 

3.00 

0.44 

1.92 

16.29 

100.00 

ORGN 

23.93 

3.47 

6.89 

13.47 

6.67 

3.87 

31.97 

3.27 

0.08 

6.39 

0.00 

100.00 

TN 

12.20 

6.39 

12.78 

16.81 

7.17 

7.70 

8.12 

3.06 

0.35 

16.52 

8.90 

100.00 

P04 

2.00 

7.13 

14.06 

9.36 

3.39 

8.79 

4.19 

1.67 

0.08 

49.33 

0.00 

100.00 

ORGP 

4.64 

7.78 

14.84 

20.38 

7.65 

8.46 

12.84 

3.75 

0.46 

19.20 

0.00 

100.00 

TP 

2.68 

6.86 

13.38 

12.13 

4.48 

8.11 

13.14 

2.20 

0.19 

36.84 

0.00 

100.00 

BOO 

28.48 

10.97 

15.21 

8.53 

6.55 

9.79 

9.06 

3.21 

0.00 

8.22 

0.00 

100.00 

SED 

16.98 

13.10 

20.98 

24.06 

6.56 

18.31 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

Segment 

200 

NH3 

2.39 

4.27 

4.95 

6.72 

0.90 

2.67 

4.86 

2.25 

0.53 

70.45 

0.00 

100.00 

N03 

15.62 

9.94 

12.14 

22.10 

6.92 

9.71 

0.43 

2.43 

0.60 

0.83 

19.29 

100.00 

ORGN 

27.83 

4.89 

4.97 

17.64 

4.81 

1.86 

32.32 

3.32 

0.14 

2.19 

0.00 

100.00 

TN 

15.97 

7.76 

9.19 

18.25 

5.30 

6.60 

7.22 

2.64 

0.49 

15.42 

11.17 

100.00 

P04 

0.87 

12.39 

15.12 

11.14 

2.98 

9.80 

5.25 

2.10 

0.17 

40.19 

0.00 

100.00 

ORGP 

5.83 

11.86 

11.54 

28.49 

5.95 

4.36 

14.00 

4.11 

0.87 

13.01 

0.00 

100.00 

TP 

2.45 

11.25 

12.76 

16.04 

3.73 

7.23 

15.51 

2.61 

0.39 

28.03 

0.00 

100.00 

BOD 

32.45 

13.77 

9.29 

11.16 

4.63 

7.20 

8.98 

3.19 

0.00 

9.38 

0.00 

100.00 

SED 

18.82 

18.84 

15.39 

30.70 

6.98 

9.25 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 
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Consequently,  we  retained  the  PQUAL  approach  used  in  the  previous  CBP  Watershed  Model 
efforts.  Atmospheric  deposition  of  N  was  added  to  the  impervious  urban  segment,  and  the 
parameters  were  adjusted  to  avoid  double-counting  of  the  atmospheric  components  in  the 
surface  runoff.  An  analogous  approach  could  not  be  implemented  for  the  urban  pervious 
segment  because  a  ’potency  factor’  approach  is  used  for  the  surface  runoff  loading  simulation, 
which  does  not  follow  a  nutrient  balance  approach  (see  Conclusions  and  Recommendations). 

e.  Estimates  of  annual  loadings  of  N  from  septic  systems  were  developed  by  CBPO  staff  (L. 
Linker,  personal  communication,  1996),  based  on  an  assessment  performed  by  Maizel  et  al. 
(1995)  using  GIS-based  procedures  as  applied  to  the  entire  Chesapeake  Bay  drainage.  The 
estimated  annual  loads  for  the  Shenandoah,  corresponding  to  10.02  lb  N/ac/yr  for  segment 
190  and  12.83  lb  N/ac/yr  for  segment  200,  were  input  as  NO3-N  as  a  constant  daily  value 
for  each  acre  of  pervious  urban  land  within  each  model  segment. 


Review  and  analysis  of  the  information  in  Table  4,  as  compared  to  the  values  from  the  earlier  study 
(Donigian  et  al.,  1995)  indicate  the  following: 

a.  The  forest  TN  and  TP  loadings  are  not  significantly  different,  primarily  because  the  earlier 
study  had  already  reduced  the  forest  TN  loadings  to  the  expected  range  of  1.0  to  2.0  lbs 
N/ac/year.  However,  the  AGCHEM  simulation  for  forests  in  this  study  does  produce  a 
somewhat  different  distribution,  with  lower  NO3-N  and  higher  Organic  N  loadings.  The  NH3- 
N  loadings  are  slightly  higher  with  the  AGCHEM  simulation,  and  the  phosphorous  simulation 
is  essentially  unchanged. 

The  BOD  loads  from  forest  have  increased  significantly  because  they  are  now  calculated  from 
the  labile  organic  N  loading  produced  by  the  enhanced  AGCHEM  module. 

b.  The  cropland  loadings  from  Hi-Till,  Lo-Till,  and  Hay  have  been  reduced  due  to  changes  in 
the  dates  and  methods  of  application  and  incorporation  of  nutrients  through  the  SPECIAL 
ACTIONS  capability.  These  changes  helped  to  reduce  some  of  the  extreme  NO3-N,  NH3-N, 
and  PO4-P  concentration  peaks  from  fertilizer  and  manure  applications,  and  thus  reduced  the 
overall  loadings  of  these  forms. 

c.  From  the  pasture  segment,  the  nutrient  simulation  with  AGCHEM  and  manure  applications 
tended  to  increase  P  loads  and  either  decrease  or  not  change  TN  loads  for  the  Shenandoah 
segments.  For  the  N  forms,  NO3-N  was  reduced  while  both  NH3-N  and  Organic  N 
increased.  While  the  relative  changes  for  NHj-N,  PO4-P,  and  Organic  P  were  significant, 
the  absolute  loadings  were  still  relatively  small,  in  the  range  of  0.2  to  0.5  Ib/ac/yr. 

d.  Loadings  from  the  Urban  Pervious  land  segment  did  not  change  significantly,  but  the  Urban 
Impervious  segment  loadings  did  increase  due  to  addition  of  the  Atmospheric  Deposition 
inputs  of  all  the  N  forms.  TN  loadings  increased  by  factors  in  the  range  of  1.5  to  2.0. 

The  distribution  of  sources  of  each  constituent  shown  in  Table  5  is  not  greatly  different  from  the 
distribution  produced  in  the  earlier  study.  Most  of  the  numbers  are  either  the  same  or  have  changed 
by  at  most  a  few  percentage  points,  which  would  not  have  a  major  impact  on  an  evaluation  of  the 
relative  contributions  of  the  different  sources.  However,  a  few  numbers  have  changed  reflecting  the 
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loading  rate  changes  noted  above  and  the  additional  sources  considered.  Some  of  these  differences  are 
as  follows: 

a.  The  forest  segments  now  represent  a  larger  relative  source  of  NOj-N  and  a  smaller  source  of 
Organic  N,  with  very  little  change  in  the  TN  percentage  contribution. 

b.  The  forest  segment  also  contributes  much  more  BOD,  due  to  the  higher  calculated  loading 
rate,  with  a  corresponding  reduction  in  the  contribution  from  other  sources  such  as  Hay,  point 
sources,  croplands,  animal  acres,  etc. 

c.  The  pasture  segment  is  now  a  larger  contributor  of  P  loads  due  to  the  higher  loading  rates 
with  the  AGCHEM  simulation,  resulting  in  a  reduction  in  the  relative  contribution  of  the  hay 
segment. 

d.  The  addition  of  the  septic  load  of  NO3-N  shows  that  it  contributes  almost  20%  of  the  load  in 
the  Shenandoah.  With  this  additional  load  and  the  reduced  forest  NO3-N  load,  the  relative 
contributions  from  other  sources  have  decreased. 

In  spite  of  these  differences,  the  key  conclusions  regarding  source  contributions  from  the  earlier 
studies  remain  unchanged.  Point  sources  remain  a  dominant  source  of  the  total  NH3-N  load  and 
represent  almost  half  of  the  total  PO4-P  load  in  the  Shenandoah.  In  terms  of  TN  and  TP,  point 
sources  represent  about  1/4  and  1/3,  respectively,  of  the  total  loads,  with  the  remaining  3/4  and  2/3 
derived  primarily  from  nonpoint  sources  and  septic  loads.  Total  agricultural  contributions  amount  to 
about  50%  of  the  TN  load  and  58%  of  the  TP  load  in  the  Shenandoah,  while  the  urban  contributions 
represent  about  10%  of  the  TN  and  6%  of  the  TP.  Although  additional  investigation  of  the  loading 
rates  is  recommended  (see  Conclusions  and  Recommendations  below),  the  results  provided  with  the 
nutrient  balance  approach  are  a  clear  improvement  from  earlier  studies  and  provide  a  sound  basis  for 
evaluating  nutrient  management  alternatives. 


Water  Quality  Calibration  Results 

Following  the  nonpoint  source  and  loading  assessment,  water  quality  calibration  involved  adjustment 
of  selected  instream  parameters  to  improve  the  overall  comparison  between  observed  and  simulated 
concentrations.  As  noted  above,  the  simulation  of  benthic  oxygen  demand  and  benthic  algae  were 
activated  in  the  earlier  effort  (Donigian  et  al.,  1995)  as  part  of  the  Phase  III  water  quality  calibration. 
These  processes  provided  additional  mechanisms  to  improve  and  better  represent  the  instream 
simulation.  The  primary  focus  of  the  calibration  was  on  maintaining  reasonable  levels  of 
phytoplankton  and  benftic  algae,  since  no  observed  data  was  available,  and  adjusting  the  related 
parameters  so  that  the  nutrients,  both  organic  and  inorganic,  followed  the  observed  variation  with  the 
revised  loadings  from  this  study.  Review  of  the  water  quality  parameters  and  the  algal  simulation  by 
the  CBPO  indicated  a  few  potential  adjustments  to  the  instream  parameters  (G.  Shenk,  personal 
communication,  1996),  and  refinements  to  the  algal  simulation  to  increase  benthic  algae 
concentrations  and  decrease  the  phytoplankton  levels  (L,  Linker,  personal  communication,  1996). 
Through  ongoing  coordination  with  CBPO  staff,  these  efforts  provided  input  to  our  calibration  work. 

Figures  6  through  13  (end  of  this  chapter)  show  the  model  results  comparisons,  in  the  form  of  eight- 
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year  timeseries  plots  of  daily  simulated  concentrations  and  the  sporadic  observed  data,  for  sediment 
(TSS),  nitrogen  forms  (NO3,  NH3,  Organic  N,  and  Total  N4),  and  phosphorous  forms  (PO4,  Organic 
P,  and  Total  P)  for  the  Shenandoah  River  at  Millville,  WV  (model  segment  200).  The  complete 
results  of  all  simulated  constituents  for  the  Millville  site  and  for  the  Shenandoah  River  at  Front  Royal, 
VA  (model  segment  190)  are  included  in  the  Appendix.  The  extent  and  volume  of  the  results 
precludes  duplication  as  part  of  this  discussion,  so  the  reader  should  refer  to  the  results  in  the 
Appendix  to  better  comprehend  the  results  discussion  below. 

The  water  quality  results  from  this  effort  are  not  greatly  different  from  those  produced  in  the  earlier 
study,  in  spite  of  all  the  changes  implemented,  including  more  extensive  application  of  AGCHEM, 
addition  of  septic  loads,  refinement  of  loading  rates,  and  additional  calibration.  However,  this  simply 
indicates  that  the  prior  calibration  was  a  good  representation  of  the  observed  data,  and  that  the 
changes  implemented  to  better  define  the  load  sources  were  incorporated  while  maintaining  the 
accuracy  of  the  overall  simulation.  The  real  benefits  of  the  current  refinement  phase  of  the  CBP 
Watershed  Model  are  realized  from  the  extension  of  the  nutrient  balance  approach  to  all  major  land 
use  (except  urban  in  our  simulations),  and  the  utility  of  this  approach  for  nutrient  management. 

Below  we  discuss  some  of  the  areas  where  the  current  results  differ  from  those  presented  earlier  in 
Donigian  et  al  (1995),  and  identify  some  problems  that  still  remain  for  selected  constituents,  where 
further  Tine  tuning’  of  the  calibration  is  recommended.  The  specific  conclusions  and  model  results 
discussed  by  Donigian  et  al  (1995)  will  not  be  repeated  here,  but  they  are  still  relevant  for  many  of 
the  water  quality  constituents  whose  simulation  did  not  change  significantly.  Some  of  the  key 
differences  in  the  simulation  results  and  issues  identified  in  this  work  are  as  follows: 

a.  The  earlier  sediment  calibration  was  generally  high  with  simulated  concentrations  higher  than 
the  sporadic  observed  values;  thus  we  attempted  to  reduce  the  concentrations  through  limited 
reductions  in  the  sediment  loading  rates.  As  shown  in  Figure  6,  this  effort  was  not  entirely 
successful  as  the  peaks  are  still  higher  than  all  of  the  observed  values  except  for  the  storm  of 
November  1985  whose  flow  was  undersimulated.  Further  reductions  in  the  loading  rates  does 
not  appear  to  be  realistic.  We  feel  that  additional  improvement  will  require  finer 
segmentation  of  the  model  to  include  shorter,  more  representative  stream  segments  to  better 
represent  localized  scour  and  deposition  processes.  Currently,  the  U.S.  Geological  Survey 
under  EPA  funding  is  developing  a  finer  segmentation  of  the  Potomac  basin  with  about  an 
order-of-magnitude  increase  in  the  number  of  land  and  stream  segments  (A.  Lumb,  personal 
communication,  1996).  Although  limited  to  the  hydrologic  and  flow  simulation,  this  more 
detailed  Potomac  model  should  also  provide  a  more  realistic  basis  for  sediment  and  water 
quality  simulation. 

b.  The  NO3-N  simulation  shown  in  Figure  7  is  acceptable  but  not  as  seasonally  correct  as  the 
earlier  results  reported  by  Donigian  et  al  (1995).  We  suspect  that  the  differences  primarily  in 
late  summer,  fall,  and  winter  for  selected  years  are  due  mostly  to  the  newly-added  septic 
loads  and  possibly  inaccuracies  in  the  seasonal  loadings  from  the  forest  and  pasture  segments, 
i.e.  the  added  AGCHEM  segments.  The  dotted  line  in  Figure  7  (  and  in  Figure  10  for  Total 
N)  shows  the  simulation  without  the  NO3-N  load  from  septic  systems  appears  to  be  improved 
in  about  five  of  the  eight  years  simulated.  As  opposed  to  being  a  constant  load,  the  septic 
load  should  probably  be  ’hydrologically’  driven  by  baseflow  so  that  the  loading  to  the  stream 
will  depend  on  the  flow  regime. 
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Both  the  NO3-N  and  Total  N  simulations  include  occasional  peaks  that  are  well  above  the 
observed  data.  This  is  typical  of  most  of  the  nutrient  simulation  results,  and  may  also  be  due 
to  the  relatively  large  model  segments  and  long  stream  reaches  used  in  the  Watershed  Model. 
However,  the  observed  data  is  relatively  infrequent  (about  every  two  months)  and  one  sample 
point  exceeding  9.0  mg/1  NO3-N  in  the  summer  of  1986  indicates  that  much  higher 
concentrations  are  possible.  Adjustments  in  the  timing  and  placement  of  nutrient  ^plications 
have  helped  to  lower  some  peaks,  but  further  investigation  and  finer  segmentation  should  be 
investigated. 

c.  Both  the  NH3-N  and  PO4-P  simulations  (Figures  8  and  1 1)  have  improved  with  reduced 
contributions  from  cropland  (along  with  the  reduced  sediment  loads),  reduced  peaks  due  to  the 
nutrient  application  adjustments,  and  increased  algal  uptake  due  to  higher  benthic  algae  levels 
than  in  the  previous  efforts.  Both  of  these  constituents  are  highly  sensitive  to  the  algal 
simulation  phytoplankton  and  benthic  algae),  for  which  no  data  was  available  for  any  of  the 
test  model  segments.  In  spite  of  the  refinements  made  in  this  effort,  algal  data  remains  as  a 
major  data  need  at  this  time;  additional  instream  benthic  and  algal  data  would  help  to  assess 
the  validity  of  the  phytoplankton  (and  benthic  algae)  simulation  and  its  impact  on  the 
inorganic  nutrient  concentrations. 

While  the  NH3-N  peaks  remain  higher  than  most  of  the  observations,  the  PO4-P  peaks  are 
much  closer  than  in  previous  efforts.  However,  base  concentrations  appear  to  be  too  low  and 
could  probably  be  improved  with  additional  calibration.  Both  of  these  constituents  would 
benefit  by  a  finer  model  segmentation  since  they  are  impacted  by  the  accuracy  of  the  sediment 
simulation. 

As  noted  earlier,  for  the  Shenandoah  the  greatest  fraction  of  the  total  loads  for  both  NH3  and 
PO4  are  derived  from  point  sources.  Thus,  although  the  concentration  peaks  in  Figures  8  and 
1 1  are  likely  due  to  nonpoint  sources,  the  point  sources  are  major  contributors  to  the  total 
load. 


c.  Simulation  of  both  Organic  N  (Figure  9)  and  Organic  P  (Figure  12)  has  improved  in  this 
effort,  and  since  Organic  P  is  the  major  component  of  Total  P  (Figure  13)  similar 
improvement  is  shown  there.  Except  for  the  November  1985  storm.  Organic  N  is  generally 
well  represented,  although  many  of  the  peaks  are  still  somewhat  high.  For  Organic  P,  the 
concentrations  are  closer  to  observed  than  in  the  previous  effort,  but  peak  concentrations  are 
also  still  high.  The  improvements  result  from  reduced  sediment  loads  (and  associated 
organics),  adjustments  to  the  N/P  ratios,  changes  to  manure  applications,  and  reductions  in 
phytoplankton  levels  (which  are  included  in  the  organic  state  variables). 

With  the  detailed  AGCHEM  simulation  of  forests,  the  forest  segment  is  now  a  major 
contributor  of  organic  N,  both  labile  and  refractory  in  both  dissolved  and  particulate  forms, 
along  with  the  animal  acres  segment.  Also,  for  the  forest  segment,  BOD  and  Organic  P 
loads  are  calculated  from  the  organic  N  components;  thus  BOD  loads  from  forest  are  also  a 
major  component.  As  shown  in  Table  5,  forests  contribute  about  25%  of  the  Organic  N  load 
and  30%  of  the  BOD  load,  while  the  animal  acres  segment  contributes  more  than  30%  of  the 
Organic  N  to  the  stream.  Both  of  these  components  need  further  investigation  as  there  is  little 
data  to  confirm  this  level  of  contribution. 
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A  bed  scour  component  for  Organic  N  was  approximated  by  CBPO,  based  on  sediment  scour 
in  the  upstream  reach  (D.  Benelmouffok,  personal  communication,  1996),  and  included  in 
their  current  modeling  efforts,  but  it  was  not  included  in  our  model  runs  due  to  a  mismatch  in 
project  schedules.  Including  bed  scour  for  organics,  in  a  manner  similar  to  current 
procedures  for  NHj  and  PO4,  is  a  relatively  minor  and  straighforward  code  enhancement  that 
should  be  implemented  to  allow  better  representation  of  major  storm  contributions  of  bed 
organics. 

Although  the  forest  simulation  provides  both  labile  and  refractory  organic  N  components 
(dissolved  and  particulate),  the  other  AGCHEM  segments  (Hi-Till,  Lo-Till,  Hay,  Pasture)  are 
restricted  to  just  the  refractory  particulate  organic  N  and  P  eroded  from  the  land  surface.  A 
more  consistent  approach  to  representing  organic  loadings  across  all  land  use  categories  is 
needed  and  should  be  investigated  jointly  with  the  algal  simulation. 


CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  results  of  our  investigation  as  presented  in  this  study,  and  our  interactions  with  CBPO 
staff  in  their  ongoing  Watershed  Model  application  as  part  of  the  Phase  IV  effort,  we  submit  the 
following  primary  conclusions  and  recommendations: 

1.  Finer  segmentation  for  all  stream  reaches  should  be  pursued  as  a  major  component  of  future 
Watershed  Model  enhancements.  Increasing  the  spatial  detail  of  the  model  by  about  a  factor 
of  lOx,  along  with  appropriate  detail  in  the  precipitation  and  land  use  inputs,  will  help  to 
improve  all  the  process  simulations,  with  specific  benefits  for  the  sediment  and  associated 
constituents,  and  benthic  processes.  Current  efforts  by  the  U.S.G.S  (restricted  to  the 
hydrology  and  flow  simulation)  should  provide  a  sound  basis  for  increasing  the  segmentation 
detail.  Similar  efforts  in  the  Patuxent  River  Basin  (AQUA  TERRA  Consultants,  1994)  may 
provide  an  indication  of  the  potential  water  quality  simulation  benefits  and  implications  of  the 
additional  segmentation  detail. 

2.  A  more  consistent  approach  for  both  BOD  and  organics  loading  needs  to  be  developed  and 
applied  consistently  across  all  land  use  categories.  Currently  the  forest  simulation 
enhancements  provide  loadings  of  both  labile  and  refractory  organic  N  components  (dissolved 
and  particulate),  while  the  other  AGCHEM  segments  (Hi-Till,  Lx)-Till,  Hay,  Pasture)  are 
restricted  to  just  the  refractory  particulate  organic  N  and  P  eroded  from  the  land  surface.  The 
forest  organic  N  capabilities  can,  and  probably  should,  be  applied  to  all  land  segments  to 
implement  this  consistent  loading  representation.  Further  investigation  of  partitioning  and 
transformation  parameters  for  the  organics  will  be  needed,  along  with  consideration  of 
extending  the  forest  N  simulation  approach  to  include  phosphorous.  In  addition,  as  noted 
above,  including  a  simple  bed  scour  algorithm  for  organics  is  a  relatively  minor  and 
straighforward  code  enhancement  that  should  be  implemented. 

3.  The  current  representation  of  septic  system  loads  needs  to  be  re-evaluated.  The  current  use  of 
a  constant  load  has  reduced  the  seasonality  of  the  NO3  simulation  shown  in  both  the  observed 
data  and  the  previous  model  simulations.  With  the  division  of  the  urban  land  use  category, 
one  alternative  approach  would  be  to  include  a  rural  residential  category  with  the  baseflow 
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concentrations  calibrated  to  match  the  annual  loading  rates  estimated  for  the  septic  systeih 
loads.  Whatever  approach  is  selected  it  should  allow  the  septic  loadings  to  be  ’hydrologically- 
driven’  so  that  the  seasonality  of  the  hydrologic  regime  and  loadings  is  represented. 

4.  The  algal  simulation,  both  phytoplankton  and  benthic  algae,  need  more  data,  investigation, 
and  evaluation.  In  this  effort,  the  benthic  algal  levels  were  increased  dramatically,  by  factors 
of  10  to  20,  based  on  very  limited  data  from  widely  scattered  sites  outside  the  Chesapeake 
Bay  drainage.  Phytoplankon  levels  showed  an  associated  decrease  from  earlier  modeling 
efforts  due  to  competition  for  the  available  nutrients.  It  seems  appropriate  that  benthic 
processes  have  a  greater  impact  than  represented  in  previous  modeling  efforts,  but  additional 
literature  data  should  be  identified  and-  actual  site-specific  data  within  the  Chesapeake  Bay 
watersheds  collected  to  confirm  the  general  magnitude  of  both  the  benthic  and  phytoplankton 
levels  represented  in  the  model.  The  algal  simulation  has  such  a  critical  impact  on  inorganic 
nutrients  levels  that  major  improvements  in  their  modeling  will  depend  on  establishing 
realistic  levels  for  the  algal  populations. 

5.  The  urban  land  use  should  be  divided  into  separate  urban  categories  —  residential, 
commercial,  industrial  -  each  with  pervious  and  impervious  fractions,  as  a  prelude  for 
nutrient  mass  balance  and  AGCHEM-type  model  application.  The  CBPO  has  pursued. an 
AGCHEM  approach  for  the  aggregated  urban  pervious  segment  in  their  current  modeling 
efforts  (L.  Linker,  personal  communication,  1996);  however,  we  feel  this  approach  requires  a 
better  definition  of  specific  urban  activities  in  order  to  develop  reasonable  nutrient  balances. 
With  a  division  into  separate  urban  categories,  nutrient  balances  would  be  more 
representative. 

6.  The  CBPO  should  explore  the  option  of  eliminating  the  current  ’composite  crop’ ' 
representation  in  the  model,  developed  as  part  of  the  Phase  II  enhancements  in  1991,  and 
move  to  simulating  each  major  crop  individually  in  each  cropland  category.  With  the 
widespread  availability  of  Pentium  processors  running  at  up  to  200  megahertz  with  gigabytes 
of  disk  storage,  many  of  the  run  time  restrictions  that  required  this  simplification  have 
disappeared.  Such  an  approach  would  allow  more  accurate  representation  of  agricultural 
practices  and  the  resulting  nutrient  balances. 

7.  In  conjunction  with  representing  each  major  crop,  the  nutrient  application  rates,  timing, 
procedures,  and  composition  distribution  (both  fertilizers  and  manure)  should  be  closely 
reviewed  and  revised  as  needed.  Experience  with  both  the  Watershed  Model  and  the  detailed 
Patuxent  Model  has  confirmed  the  critical  importance  of  the  assumptions  underlying  the 
nutrient  applications  in  the  model.  Issues  such  as  double-cropping,  application  methods, 
inorganic/organic  splits  (e.g.  VOJOtg  P  splits  for  manure  and  sludge),  etc.  need  to  be 
investigated  to  improve  the  nutrient  mass  balances  and  simulations. 

8.  The  forest  N  simulation  approach  should  be  extended  to  include  the  P  cycle,  so  that  both  N 
and  P  mass  balances  can  be  implemented  for  all  land  segments.  Just  as  the  previous 
AGCHEM  module  provided  a  valid  framework  for  the  forest  N  enhancements,  the  P  cycle 
processes  currently  in  AGCHEM  can  be  readily  adapted  for  forested  conditions.  Field  site 
testing  on  small  forested  watersheds  would  be  needed  to  fiilly  evaluate  the  code  enhancements. 
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Appendix  A 

Shenandoah  Model  Segments  Results 

Simulated  and  Observed  Flow  at  Reach  190, 1984-87 

Simulated  and  Observed  Flow  at  Reach  190, 1988-91 

Frequency  Analysis  of  Flow  at  Reach  190 

Comparison  of  Annual  Total  Observed  and  Simulated  Flow  Volume 

Simulated  Se<fiment(TSS)  Concentration  at  Reach  190 

Simulated  and  Observed  Water  Temperature  at  Reach  190 

Simulated  and  Observed  DO  Concentration  at  Reach  190 

Simulated  and  Observed  Nitrate-N  Concentration  at  Reach  190 

Simulated  Ammonia-N  Concentration  at  Reach  190 

Simulated  Organic-N  Concentration  at  Reach  190 

Simulated  Total-N  Concentration  at  Reach  190 

Simulated  and  Observed  PO4-P  Concentration  at  Reach  190 

Simulated  Organic-P  Concentration  at  Reach  190 

Simulated  Total-P  Concentration  at  Reach  190 

Simulated  TOC  Concentration  at  Reach  190 

Simulated  Chlorophyll  A  Concentration  at  Reach  190 

Simulated  BOD  Concentration  at  Reach  190 

Simulated  Benthic  Algae  Concentration  at  Reach  190 

Simulated  and  Observed  Flow  at  Reach  200, 1984-87 
Simulated  and  Observed  Flow  at  Reach  200, 1988-91 
Frequency  Analysis  of  Flow  at  Reach  200 

Simulated  and  Observed  Sediment(TSS)  Concentration  at  Reach  200 
Simulated  and  Observed  Water  Temperature  at  Reach  200 
Simulated  and  Observed  DO  Concentration  at  Reach  200 
Simidated  and  Observed  Nitrate-N  Concentration  at  Reach  200 
Simulated  and  Observed  Ammonia-N  Concentration  at  Reach  200 
Simulated  and  Observed  Organic-N  Concentration  at  Reach  200 
Simulated  and  Observed  Total-N  Concentration  at  Reach  200 
Simulated  and  Observed  PO4-P  Concentration  at  Reach  200 
Simulated  and  Observed  Organic-P  Concentration  at  Reach  200 
Simulated  and  Observed  Total-P  Concentration  at  Reach  200 
Simulated  TOC  Concentration  at  Reach  200 
Simulated  Chlorophyll  A  Concaitration  at  Reach  200 
Simulated  BOD  Concentration  at  Reach  200 


Appendix  A  Shenandoah  Model  Segments  Results 


A1 


Simulated  Benthic  Algae  Concentration  at  Reach  200 
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AGCHEM  Summary  for  Shenandoah  Basin  (Low-Till),  PERLND  193 
AGCHEM  Summary  for  Shenandoah  Basin  (PASTURE),  PERLND  194 
AGCHEM  Summary  for  Shenandoah  Basin  (HAY),  PERLND  196 
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Per  Acre  Load  Contributed  from  Each  Land  Use  in  Shenandoah  Basin  (Ib/ac) 
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1984  1985  1986  1987 

SIMULATED  AND  OBSERVED  FLOW  AT  REACH  190 
SF  SHENANDOAH  RIVER  AT  FRONT  ROYAL,  VA 


D  AND  OBSERV 
NDOAH  RIVER 


CHESAPEAKE  BAY  WATERSHED  HYDROLOGIC  CALIBRA  TION 
COMPARISON  OF  ANNUAL  TOTAL  OBSERVED  vs  SIMULATED  FLOW 

PHASE  IV 


SF  SHENANDOAH  RIVER  AT  FRONT  ROYAL,  VA  (SEGMENT  190) 


YEAR 

OBSERVED* 

SIMULATED 

FLOW  (in) 

FLOW  (in) 

1984 

18.60 

20.05 

198.5 

16.20 

17.10 

1986 

7..34 

7.24 

1987 

16.10 

14.73 

1988 

7,71 

7.06 

1989 

15.08 

16..35 

1990 

l.^JiS 

15.41 

1991 

11.67 

11.63 

MEAN 

13.29 

13.70 

*  -  Observed  Flow  at  SF  Shenandoah  River  at  Front  Royal.  VA 
**  -  Simulated  Outflow  from  RCH  190 


SHENANDOAH  RIVER  AT  MILLVILLE,  WV  (SEGMENT  200) 


YEAR 

OBSERVED* 

SIMULATED 

FLOW  (in) 

FLOW  (in) 

1984 

17,82 

18.27 

1985 

12.99 

14..35 

1986 

6.92 

6.88 

1987 

13.95 

12.09 

1988 

7.72 

6.72 

1989 

12.02 

13.15 

1990 

11.58 

13.78 

1991 

10.50 

9.99 

MEAN 

11.69 

11.90 

*  -  Observed  Flow  at  Shenandoah  River  at  Millville,  WV 
**  -  Simulated  Outflow  from  RCH  200 
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PER  ACRE  LOAD  CONTRIBUTED  FROM  EACH  LAND  USE  IN  SHENANDOAH  BASIN  (LB/AC);  SEDIMENT  IN  TONS/AC 


Segment  190 


•Pervious- 

-X - Impervious-- 

•>  Total 

FOR 

HTC 

LTC 

PAS 

URB 

NAY 

ANML 

RES 

Load 

NH3 

1984 

0.132 

4.741 

5.016 

0.896 

0.521 

0.791 

268.145 

2.320 

2.204 

1985 

0.086 

1.756 

1.254 

0.561 

0.372 

0.411 

264.176 

2.090 

1.607 

1986 

0.052 

1.907 

1.462 

0.241 

0.129 

0.209 

145.689 

1.800 

1.454 

1987 

0.088 

2.964 

2.354 

0.591 

0.394 

0.430 

255.534 

2.260 

1.837 

1988 

0.040 

1.649 

1.455 

0.221 

0.114 

0.172 

156.103 

1.880 

1.590 

1989 

0.075 

3.680 

2.824 

0.740 

0.329 

1.050 

256.157 

2.700 

1.713 

1990 

0.078 

3.054 

2.649 

0.607 

0.372 

0.464 

235.429 

2.220 

1.575 

1991 

0.077 

2.053 

1.566 

0.497 

0.231 

0.435 

190.798 

1.890 

1.374 

MEAN 

0.078 

2.726 

2.323 

0.544 

0.308 

0.495 

221.504 

2.145 

1.669 

NH3 

1984 

1.772 

17.710 

15.390 

4.777 

7.940 

5.359 

67.036 

6.670 

5.316 

1985 

0.953 

17.310 

15.230 

4.103 

6.690 

3.620 

66.044 

5.960 

4.420 

1986 

0.715 

11.570 

8.992 

3.150 

3.490 

2.492 

36.422 

5.740 

3.247 

1987 

1.282 

16.800 

13.200 

4.969 

6.240 

4.470 

63.883 

6.560 

4.743 

1988 

0.437 

14.860 

12.450 

3.206 

3.430 

2.855 

39.026 

5.590 

3.412 

1989 

0.736 

24.080 

21.870 

5.559 

6.940 

5.321 

64.039 

7.340 

5.362 

1990 

0.524 

21.170 

19.790 

4.778 

6.490 

4.299 

58.857 

6.570 

4.741 

1991 

0.622 

14.420 

12.210 

4.180 

5.080 

3.556 

47.700 

5.900 

3.909 

MEAN 

0.880 

17.240 

14.892 

4.340 

5.788 

3.997 

55.376 

6.291 

4.394 

ORGN 

1984 

1.492 

5.649 

5.048 

2.143 

3.569 

1.279  2011.084 

3.395 

3.035 

1985 

1.218 

3.961 

3.657 

1.394 

2.445 

0.998 

1981.320 

3.061 

2.478 

1986 

0.608 

0.981 

0.757 

0.652 

0.675 

0.167 

1092.665 

3.235 

1.196 

1987 

1.129 

5.327 

4.470 

1.476 

2.675 

1.037 

1916.503 

3.335 

2.531 

1988 

0.609 

1.454 

1.161 

0.482 

0.549 

0.117  1170.770 

3.150 

1.235 

1989 

1.100 

2.718 

2.147 

1,620 

2.033 

0.557  1921.181 

3,558 

2.279 

1990 

1.174 

3.869 

3.423 

1.414 

2.452 

0.947 

1765.716 

3.376 

2.372 

1991 

0.828 

2.632 

2.069 

1.244 

1.380 

0.365 

1430.987 

3.231 

1.778 

MEAN 

1.020 

3.324 

2.842 

1.303 

1.972 

0.683 

1661.278 

3.293 

2.113 

TN 

1984 

3.395 

28.100 

25.460 

7.816 

12.030 

7.429 

1944.048 

12.385 

10.392 

1985 

2.257 

23.020 

20.140 

6.058 

9.507 

5.028 

1915,276 

11.111 

8.343 

1986 

1.375 

14.460 

11.210 

4.044 

4.294 

2.867 

1056.243 

10.775 

5.808 

1987 

2.499 

25.090 

20.020 

7.035 

9.309 

5.937  1852.619 

12.155 

8.955 

1988 

1.085 

17.970 

15.070 

3.909 

4.093 

3.144 

1131.745 

10.620 

6.141 

1989 

1.910 

30.480 

26.840 

7.919 

9.302 

6.928 

1857.141 

13.598 

9.198 

1990 

1.776 

28.090 

25.860 

6.799 

9.314 

5.709 

1706.858 

12.166 

8,545 

1991 

1.527 

19.100 

15.840 

5.921 

6.691 

4.355 

1383.287 

11.021 

6.944 

MEAN 

1.978 

23.289 

20.055 

6.188 

8.068 

5.175 

1605.902 

11.729 

8.041 

P04 

1984 

0.036 

1.559 

1.906 

0.249 

0.472 

0.548 

67.036 

0.445 

0.687 

1985 

0.021 

1.846 

1.291 

0.266 

0.321 

0.347 

66.044 

0.396 

0.593 

1986 

0.006 

1.012 

0.788 

0.131 

0.079 

0.158 

36.422 

0.418 

0.430 

1987 

0.033 

2.015 

1.556 

0.280 

0.353 

0.469 

63.883 

0.434 

0.654 

1988 

0.008 

1.270 

1.150 

0.068 

0.061 

0.086 

39.026 

0.406 

0.427 

1989 

0.026 

2.761 

2.333 

0.284 

0.261 

0.646 

64.039 

0.469 

0.568 

1990 

0.026 

2.230 

1.811 

0.321 

0.321 

0.420 

58.857 

0.442 

0.513 

1991 

0.018 

1.189 

0.969 

0.244 

0.173 

0.488 

47.700 

0.419 

0.429 

MEAN 

0.022 

1.735 

1.475 

0.230 

0.255 

0.395 

55.376 

0.429 

0.538 

ORGP 

1984 

0.036 

1.571 

1.348 

0.405 

0.510 

0.343 

100.554 

0.485 

0.401 

1985 

0.028 

1.109 

0.985 

0.262 

0.349 

0.275 

99.066 

0.437 

0.311 

1986 

0.014 

0.273 

0.203 

0.123 

0.096 

0.046 

54.633 

0.462 

0.146 

1987 

0.029 

1.486 

1.198 

0.279 

0.382 

0.281 

95.825 

0.476 

0.340 

1988 

0.015 

0.405 

0.311 

0.088 

0.078 

0.031 

58.539 

0,450 

0.158 

1989 

0.026 

0.759 

0.577 

0.305 

0.290 

0.153 

96.059 

0.508 

0.253 

1990 

0.030 

1.080 

0.919 

0.266 

0.350 

0.261 

88.286 

0.482 

0.285 

1991 

0.020 

0.734 

0.554 

0.236 

0.197 

0.099 

71.549 

0.462 

0.211 

MEAN 

0.025 

0.927 

0.762 

0.245 

0.282 

0.186 

83.064 

0.470 

0.263 
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Segment  190 


•Pervious-' 

•X - Impervious- -> 

Total 

FOR 

HTC 

LTC 

PAS 

URB 

HAY 

ANML 

RES 

Load 

TP 

1984 

0.071 

3.130 

3.254 

0.654 

0.982 

0.891 

335.181 

0.930 

1.156 

1985 

0.049 

2.956 

2.276 

0.528 

0.670 

0.622 

330.220 

0.833 

0.971 

1986 

0.020 

1.285 

0.990 

0.254 

0.176 

0.204 

182.111 

0.880 

0.613 

1987 

0.061 

3.501 

2.755 

0.559 

0.735 

0.750 

319.417 

0.910 

1.058 

1988 

0.023 

1.675 

1.460 

0.156 

0.139 

0.117 

195.128 

0.856 

0.625 

1989 

0.052 

3.520 

2.910 

0.588 

0.551 

0.799 

320.197 

0.977 

0.887 

1990 

0.056 

3.310 

2.730 

0.587 

0.671 

0.680 

294.286 

0.924 

0.857 

1991 

0.039 

1.922 

1.523 

0.480 

0.370 

0.587 

238.498 

0.881 

0.689 

MEAN 

0.046 

2.662 

2.237 

0.476 

0.537 

0.581 

276.880 

0.899 

0.857 

BOD 

1984 

14.795 

123.300 

74.700 

11.010 

28.860 

25.200  4692.528 

27.450 

25.313 

1985 

12.537 

116.100 

72.300 

7.350 

19.770 

20.490  4623.079 

24.750 

21.486 

1986 

6.305 

20.520 

11.730 

3.450 

5.460 

3.810  2549.551 

26.160 

8.551 

1987 

10.579 

141.300 

80.700 

7.650 

21.630 

21.000  4471.839 

26.970 

22.137 

1988 

5.959 

32.100 

23.490 

2.646 

4.440 

3.120  2731.798 

25.470 

9.315 

1989 

10.832 

75.000 

45.900 

8.460 

i6.ao 

12.690  4482.755 

28.770 

17.875 

1990 

10.940 

113.700 

68.400 

7.410 

19.830 

19.350  4120.003 

27.300 

20.567 

1991 

8.009 

69.900 

36.000 

6.390 

11.160 

8.280  3338.969 

26.130 

13.978 

MEAN 

9.995 

86.490 

51.653 

6.796 

15.949 

14.243  3876.315 

26.625 

17.403 

SED 

1984 

0.097 

1.680 

1.210 

0.305 

0.291 

0.469 

0.000 

0.000 

0.291 

1985 

0.058 

1.160 

0.864 

0.193 

0.189 

0.366 

0.000 

0.000 

0.198 

1986 

0.018 

0.287 

0.178 

0.089 

0.043 

0.062 

0.000 

0.000 

0.054 

1987 

0.078 

1.570 

1.060 

0.207 

0.215 

0.377 

0.000 

0.000 

0.233 

1988 

0.025 

0.429 

0.275 

0.064 

0.032 

0.046 

0.000 

0.000 

0.058 

1989 

0.056 

0.806 

0.517 

0.226 

0.149 

0.213 

0.000 

0.000 

0.157 

1990 

0.081 

1.130 

0.811 

0.197 

0.193 

0.346 

0.000 

0.000 

0.205 

1991 

0.041 

0.772 

0.489 

0.173 

0.100 

0.142 

0.000 

0.000 

0.124 

MEAN 

0.057 

0.979 

0.676 

0.182 

0.152 

0.253 

0.000 

0.000 

0.165 
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Segment  200 


<• 

FOR 

HTC 

LTC 

•Pervious-- 

PAS 

URB 

HAY 

•X - Impervious--> 

ANHL  RES 

Total 

Load 

NH3 

1984 

0.106 

3.142 

2.695 

0.594 

0.366 

0.584 

234.778 

2.170 

1.636 

1985 

0.061 

1.334 

0.895 

0.351 

0.214 

0.235 

211.834 

1.910 

1.405 

1986 

0.041 

2.172 

1.583 

0.204 

0.133 

0.238 

129.661 

1.680 

1.320 

1987 

0.059 

2.178 

1.632 

0.342 

0.214 

0.266 

214.257 

2.020 

1.418 

1988 

0.033 

1.715 

1.574 

0.219 

0.121 

0.181 

139.478 

1.700 

1.366 

1989 

0.048 

2.165 

1.810 

0.409 

0.232 

0.422 

194.622 

2.390 

0.970 

1990 

0.057 

2.166 

2.032 

0.446 

0.294 

0.340 

223.523 

2.140 

0.992 

1991 

0.046 

1.183 

0.839 

0.259 

0.120 

0.274 

134.624 

1.580 

0.847 

MEAN 

0.057 

2.007 

1.632 

0.353 

0.212 

0.317 

185.347 

1.949 

1.244 

N03 

1984 

1.584 

15.050 

12.790 

4.141 

6.790 

4.868 

58.694 

6.380 

4.411 

1985 

0.804 

11.550 

10.400 

2.935 

5.080 

2.595 

52.959 

5.700 

3.188 

1986 

0.987 

9.438 

7.347 

2.784 

3.390 

2.377 

32.415 

5.490 

2.948 

1987 

2.940 

13.040 

9.713 

3.499 

4.440 

3.365 

53.564 

6.140 

4.513 

1988 

0.497 

14.170 

12.130 

2.647 

3.210 

2.530 

34.870 

5.180 

2.965 

1989 

0.630 

18.930 

16.690 

3.541 

4.690 

3.789 

48.656 

6.780 

3.788 

1990 

0.418 

15.010 

15.350 

3.792 

5.590 

3.707 

55.881 

6.450 

3.618 

1991 

0.478 

8.301 

6.020 

2.908 

3.590 

2.809 

33.656 

5.390 

2.683 

MEAN 

1.042 

13.186 

11.305 

3.281 

4.598 

3.255 

46.337 

5.939 

3.514 

ORGN 


1984 

1.224 

6.451 

4.966 

2.004 

2.360 

0.744 

1760.834 

3.317 

2.337 

1985 

0.839 

2.167 

1.788 

0.983 

1.254 

0.260 

1588.758 

3.072 

1.491 

1986 

0.474 

0.723 

0.481 

0.630 

0.720 

0.045 

972.455 

3.168 

0.872 

1987 

0.593 

1.998 

1.263 

1.121 

1.317 

0.119 

1606.928 

3.276 

1.366 

1988 

0.438 

1.278 

0.548 

0.418 

0.646 

0.031 

1046.086 

3.020 

0.840 

1989 

0.600 

2.421 

1.537 

1.203 

1.451 

0.145 

1459.666 

3.473 

1.373 

1990 

1.203 

4.053 

2.969 

1.579 

1.881 

0.481 

1676.423 

3.376 

2.005 

1991 

0.563 

1.652 

1.247 

0.427 

0.582 

0.164 

1009.678 

3.187 

0.942 

MEAN 

0.742 

2.593 

1.850 

1.046 

1.276 

0.249 

1390.104 

3.236 

1.403 

TN 

1984 

2.914 

24.650 

20.450 

6.740 

9.516 

6.196 

1702.139 

11.867 

8.269 

1985 

1.704 

15.050 

13.080 

4.269 

6.548 

3.090 

1535.799 

10.682 

5.980 

1986 

1.503 

12.330 

9.410 

3.618 

4.243 

2.660 

940.040 

10.338 

5.077 

1987 

3.592 

17.210 

12.610 

4.963 

5.971 

3.750 

1553.364 

11.436 

7.192 

1988 

0.967 

17.160 

14.250 

3.285 

3.977 

2.742 

1011.217 

9.900 

5.103 

1989 

1.277 

23.520 

20.030 

5.152 

6.373 

4.356 

1411.011 

12.643 

6.035 

1990 

1.679 

21.220 

20.350 

5.817 

7.765 

4.528 

1620.542 

11.966 

6.505 

1991 

1.088 

11.140 

8.107 

3.594 

4.292 

3.247 

976.022 

10.157 

4.406 

MEAN 

1.841 

17.785 

14.786 

4.680 

6.085 

3.821 

1343.767 

11.124 

6.071 

P04 

1984 

0.010 

1.015 

1.040 

0.182 

0.308 

0.423 

58.694 

0.435 

0.326 

1985 

0.005 

1.374 

1.009 

0.149 

0.159 

0.160 

52.959 

0.397 

0.287 

1986 

0.001 

1.045 

0.752 

0.091 

0.086 

0.176 

32.415 

0.409 

0.235 

1987 

0.002 

1.464 

1.135 

0.147 

0.169 

0.285 

53.564 

0.426 

0.307 

1988 

0.001 

1.425 

1.372 

0.061 

0.077 

0.089 

34.870 

0.387 

0.265 

1989 

0.003 

2.093 

1.992 

0.130 

0.187 

0.357 

48.656 

0.457 

0.330 

1990 

0.012 

1.631 

1.393 

0.205 

0.245 

0.389 

55.881 

0.442 

0.326 

1991 

0.003 

0.717 

0.527 

0.118 

0.066 

0.273 

33.656 

0.412 

0.225 

MEAN 

0.005 

1.345 

1.153 

0.135 

0.162 

0.269 

46.337 

0.421 

0.288 

ORGP 

1984 

0.030 

1.802 

1.328 

0.377 

0.337 

0.200 

88.042 

0.474 

0.301 

1985 

0.020 

0.608 

0.481 

0.181 

0.179 

0.071 

79.438 

0.439 

0.161 

1986 

0.010 

0.201 

0.129 

0.117 

0.103 

0.012 

48.623 

0.453 

0.094 

1987 

0.013 

0.559 

0.338 

0.210 

0.188 

0.032 

80.346 

0.468 

0.154 

1988 

0.009 

0.358 

0.147 

0.075 

0.092 

0.008 

52.304 

0.431 

0.094 

1989 

0.014 

0.677 

0.412 

0.225 

0.207 

0.039 

72.983 

0.496 

0.159 

1990 

0.035 

1.138 

0.797 

0.298 

0.269 

0.131 

83.821 

0.482 

0.231 

1991 

0.013 

0.462 

0.333 

0.077 

0.083 

0.044 

50.484 

0.455 

0.103 

MEAN 

0.018 

0.726 

0.496 

0.195 

0.182 

0.067 

69.505 

0.462 

0.162 
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Segment  200 


< - - 

•Pervious-' 

■X - Impervious--> 

Total 

FOR 

HTC 

UTC 

PAS 

UR8 

HAY 

ANHL 

RES 

Load 

TP 

1984 

0.040 

2.817 

2.368 

0.558 

0.645 

0.623 

293.472 

0.909 

0.674 

1985 

0.025 

1.982 

1.490 

0.331 

0.338 

0.232 

264.793 

0.836 

0.491 

1986 

0.011 

1.246 

0.881 

0.208 

0.189 

0.189 

162.076 

0.862 

0.356 

1987 

0.016 

2.023 

1.473 

0.357 

0.357 

0.317 

267.821 

0.894 

0.505 

1988 

0.010 

1.783 

1.519 

0.137 

0.169 

0.098 

174.348 

0.818 

0.387 

1989 

0.017 

2.770 

2.404 

0.355 

0.394 

0.397 

243.278 

0.953 

0.529 

1990 

0.047 

2.769 

2.189 

0.503 

0.514 

0.520 

279.404 

0.924 

0.602 

1991 

0.016 

1.178 

0.860 

0.194 

0.149 

0.317 

168.280 

0.867 

0.356 

MEAN 

0.023 

2.071 

1.648 

0.330 

0.344 

0.336 

231.684 

0.883 

0.488 

BOO 

1984 

11.932 

144.000 

66.900 

10.410 

19.080 

19.830  4108.612 

26.820 

20.921 

1985 

8.383 

58.800 

32.100 

5.370 

10.140 

9.750  3707.101 

24.840 

12.268 

1986 

5.249 

15.150 

7.890 

3.390 

5.820 

1.779  2269.063 

25.620 

6.727 

1987 

6.283 

45.600 

20.160 

5.910 

10.650 

4.470  3749.498 

26.490 

10.566 

1988 

4.759 

26.340 

11.040 

2.379 

5.220 

1.974  2440.868 

24.420 

7.074 

1989 

6.214 

57.300 

28.140 

6.330 

11.730 

6.240  3405.888 

28.080 

11.397 

1990 

9.549 

111.000 

50.700 

8.220 

15.210 

16.620  3911.653 

27.300 

17.299 

1991 

5.747 

32.100 

15.120 

2.451 

4.710 

4.140  2355.914 

25.770 

8.106 

MEAN 

7.265 

61.286 

29.006 

5.558 

10.320 

8.100  3243.574 

26.168 

11.795 

SED 

1984 

0.082 

1.920 

1.180 

0.274 

0.280 

0.275 

0.000 

0.000 

0.247 

1985 

0.045 

0.640 

0.429 

0.131 

0.133 

0.102 

0.000 

0.000 

0.105 

1986 

0.003 

0.215 

0.116 

0.083 

0.074 

0.019 

0.000 

0.000 

0.037 

1987 

0.012 

0.598 

0.304 

0.151 

0.149 

0.051 

0.000 

0.000 

0.083 

1988 

0.001 

0.380 

0.134 

0.053 

0.065 

0.016 

0.000 

0.000 

0.034 

1989 

0.014 

0.715 

0.370 

0.163 

0.168 

0.060 

0.000 

0.000 

0.094 

1990 

0.137 

1.200 

0.701 

0.215 

0.222 

0.177 

0.000 

0.000 

0.212 

1991 

0.015 

0.489 

0.296 

0.052 

0.052 

0.063 

0.000 

0.000 

0.054 

MEAN 

0.039 

0.770 

0.441 

0.140 

0.143 

0.095 

0.000 

0.000 

0.108 
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PERCENT  OF  TOTAL  LOAD  CONTRIBUTED  FROH  EACH  LAND  USE/SOURCE  IN  SHENANDOAH  BASIN 


Segment  190 


<- - 

FOR 

HTC 

---Pervic 

LTC 

HJS - - 

PAS 

URB 

HAY 

-X - Impervious-- 

ANML  RES 

->  Atmos 
Dep 

Point 

Source 

Septic 

Load 

Total 

Load 

NH3 

1984 

2.96 

4.75 

11.66 

8.88 

1.69 

4.29 

4.95 

2.21 

0.37 

58.24 

0.00 

100.00 

1985 

2.65 

2.41 

4.00 

7.63 

1.65 

3.06 

6.68 

2.72 

0.42 

68.75 

0.00 

100.00 

1986 

1.77 

2.90 

5.15 

3.63 

0.63 

1.72 

4.08 

2.59 

0.27 

77.28 

0.00 

100.00 

1987 

2.38 

3.56 

6.57 

7.02 

1.53 

2.80 

5.66 

2.58 

0.41 

67.51 

0.00 

100.00 

1988 

1.24 

2.29 

4.69 

3.03 

0.51 

1.30 

3.99 

2.48 

0.30 

80.15 

0.00 

100.00 

1989 

2.17 

4.74 

8.44 

9.44 

1.37 

7.33 

6.08 

3.30 

0.52 

56.60 

0.00 

100.00 

1990 

2.46 

4.28 

8.62 

8.42 

1.69 

3.52 

6.08 

2.95 

0.44 

61.56 

0.00 

100.00 

1991 

2.79 

3.30 

5.84 

7.91 

1.20 

3.78 

5.65 

2.88 

0.41 

66.25 

0.00 

100.00 

MEAN 

2.33 

3.60 

7.13 

7.12 

1.32 

3.55 

5.39 

2.69 

0.39 

66.45 

0.00 

100.00 

NOS 

1984 

16.53 

7.36 

14.83 

19.64 

10.67 

12.06 

0.51 

2.63 

0.44 

1.88 

13.46 

100.00 

1985 

10.69 

8.65 

17.66 

20.28 

10.82 

9.80 

0.61 

2.82 

0.47 

2.02 

16.20 

100.00 

1986 

10.91 

7.87 

14.19 

21.19 

7.68 

9.18 

0.46 

3.70 

0.43 

2.35 

22.04 

100.00 

1987 

13.40 

7.82 

14.26 

22.88 

9.40 

11.27 

0.55 

2.90 

0.45 

1.98 

15.09 

100.00 

1988 

6.34 

9.62 

18.70 

20.53 

7.18 

10.01 

0.47 

3.43 

0.44 

2.31 

20.98 

100.00 

1989 

6.80 

9.92 

20.89 

22.64 

9.25 

11.87 

0.49 

2.87 

0.42 

1.49 

13.35 

100.00 

1990 

5.48 

9.86 

21.38 

22.02 

9.78 

10.85 

0.50 

2.90 

0.44 

1.68 

15.10 

100.00 

1991 

7.89 

8.14 

16.01 

23.36 

9.29 

10.88 

0.50 

3.16 

0.45 

2.00 

18.31 

100.00 

MEAN 

9.93 

8.66 

17.36 

21.58 

9.41 

10.88 

0.51 

3.00 

0.44 

1.92 

16.29 

100.00 

ORGN 

1984 

24.38 

4.11 

8.52 

15.42 

8.40 

5.04 

26.94 

2.34 

0.07 

4.78 

0.00 

100.00 

■1985 

24.38 

3.53 

7.56 

12.29 

7.05 

4.82 

32.51 

2.59 

0.07 

5.21 

0.00 

100.00 

1986 

25.21 

1.81 

3.24 

11.91 

4.03 

1.67 

37.13 

5.67 

0.12 

9.22 

0.00 

100.00 

1987 

22.13 

4.65 

9.05 

12.74 

7.55 

4.90 

30.78 

2.76 

0.07 

5.38 

0.00 

100.00 

1988 

24.45 

2.60 

4.82 

8.52 

3.18 

1.13 

38.54 

5.34 

0.12 

11.27 

0.00 

100.00 

1989 

23.92 

2.63 

4.83 

15.52 

6.38 

2.92 

34.27 

3.27 

0.08 

6.19 

0.00 

100.00 

1990 

24.53 

3.60 

7.39 

13.02 

7.39 

4.78 

30.27 

2.98 

0.08 

5.95 

0.00 

100.00 

1991 

23.08 

3.27 

5.96 

15.28 

5.55 

2.46 

32.72 

3.81 

0.09 

7.80 

0.00 

100.00 

MEAN 

23.93 

3.47 

6.89 

13.47 

6.67 

3.87 

31.97 

3.27 

0.08 

6.39 

0.00 

100.00 

TN 


1984 

16.20 

5.97 

12.56 

16.43 

8.27 

8.55 

7.61 

2.50 

0.32 

14.71 

6.89 

100.00 

1985 

13.41 

6.09 

12.37 

15.86 

8.14 

7.21 

9.33 

2.79 

0.35 

15.86 

8.58 

100.00 

1986 

11.74 

5.50 

9.89 

15.21 

5.28 

5.90 

7.39 

3.89 

0.33 

22.54 

12.32 

100.00 

1987 

13.84 

6.19 

11.45 

17.16 

7.43 

7.93 

8.41 

2.84 

0.34 

16.41 

7.99 

100.00 

1988 

8.76 

6.46 

12.57 

13.90 

4.76 

6.12 

7.49 

3.62 

0.34 

24.31 

11.65 

100.00 

1989 

10.30 

7.32 

14.95 

18.80 

7.23 

9.01 

8.21 

3.10 

0.36 

12.94 

7.78 

100.00 

1990 

10.31 

7.26 

15.50 

17.38 

7.79 

7.99 

8.12 

2.98 

0.35 

13.93 

8.38 

100.00 

1991 

10.90 

6.07 

11.69 

18.62 

6.89 

7.50 

8.10 

3.32 

0.36 

16.23 

10.31 

100.00 

MEAN 

12.20 

6.39 

12.78 

16.81 

7.17 

7.70 

8.12 

3.06 

0.35 

16.52 

8.90 

100.00 

P04 

1984 

2.57 

5.01 

14.22 

7.92 

4.91 

9.55 

3.97 

1.36 

0.07 

50.43 

0.00 

100.00 

1985 

1.74 

6.88 

11.16 

9.80 

3.87 

6.99 

4.53 

1.40 

0.07 

53.56 

0.00 

100.00 

1986 

0.70 

5.20 

9.38 

6.66 

1.32 

4.39 

3.44 

2.04 

0.08 

66.79 

0.00 

100.00 

1987 

2.48 

6.81 

12.20 

9.36 

3.86 

8.59 

3.97 

1.39 

0.07 

51.27 

0.00 

100.00 

1988 

0.94 

6.57 

13.80 

3.46 

1.02 

2.41 

3.72 

1.99 

0.08 

66.03 

0.00 

100.00 

1989 

2.26 

10.73 

21.04 

10.90 

3.28 

13.60 

4.58 

1.73 

0.08 

31.81 

0.00 

100.00 

1990 

2.55 

9.60 

18.10 

13.66 

4.48 

9.79 

4.67 

1.81 

0.09 

35.26 

0.00 

100.00 

1991 

2.13 

6.11 

11.56 

12.42 

2.88 

13.59 

4.52 

2.04 

0.09 

44.66 

0.00 

100.00 

MEAN 

2.00 

7.13 

14.06 

9.36 

3.39 

8.79 

4.19 

1.67 

0.08 

49.33 

0.00 

100.00 

ORGP 

1984 

4.40 

8.65 

17.23 

22.06 

9.09 

10.24 

10.20 

2.53 

0.35 

15.24 

0.00 

100.00 

1985 

4.45 

7.88 

16.25 

18.42 

8.04 

10.60 

12.97 

2.95 

0.41 

18.04 

0.00 

100.00 

1986 

4.68 

4.13 

7.11 

18.35 

4.72 

3.75 

15.21 

6.63 

0.67 

34.72 

0.00 

100.00 

1987 

4.17 

9.66 

18.08 

17.95 

8.04 

9.90 

11.47 

2.94 

0.37 

17.42 

0.00 

100.00 

1988 

4.62 

5.66 

10.08 

12.24 

3.55 

2.37 

15.06 

5.97 

0.64 

39.79 

0.00 

100.00 

1989 

5.16 

6.61 

11.66 

26.25 

8.19 

7.21 

15.41 

4.20 

0.52 

14.78 

0.00 

100.00 

1990 

5.19 

8.37 

16.53 

20.42 

8.79 

10.94 

12.60 

3.54 

0.46 

13.15 

0.00 

100.00 

1991 

4.75 

7.66 

13.44 

24.39 

6.67 

5.61 

13.77 

4.58 

0.55 

18.60 

0.00 

100.00 

MEAN 

4.64 

7.78 

14.84 

20.38 

7.65 

8.46 

12.84 

3.75 

0.46 

19.20 

0.00 

100.00 
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Segment  190 


<• 

---Pervu 

>US . 

-X - Impervious- 

--->  Atnx 

5S  Point 

Septic 

Total 

FOR 

HTC 

LTC 

PAS 

URB 

HAY 

ANML 

RES 

Dep 

Source 

Load 

Load 

TP 

1984 

3.05 

5.98 

14.43 

12.36 

6.07 

9.22 

11.79 

1.69 

0.16 

35.25 

0.00 

100.00 

1985 

2.49 

6.72 

12.02 

11.88 

4.94 

7.66 

13.83 

1.80 

0.17 

38.49 

0.00 

100.00 

1986 

1.61 

4.63 

8.28 

9.04 

2.05 

3.97 

12.08 

3.01 

0.21 

55.12 

0.00 

100.00 

1987 

2.87 

7.30 

13.34 

11.54 

4.96 

8.48 

12.27 

1.80 

0.16 

37.24 

0.00 

100.00 

1988 

1.81 

5.92 

11.98 

5.46 

1.60 

2.25 

12.70 

2.87 

0.22 

55.20 

0.00 

100.00 

1989 

2.92 

8.76 

16.81 

14.49 

4.44 

10.78 

14.68 

2.31 

0.20 

24.61 

0.00 

100.00 

1990 

3.25 

8.52 

16.32 

14.96 

5.60 

9.49 

13.95 

2.26 

0.21 

25.45 

0.00 

100.00 

1991 

2.78 

6.16 

11.32 

15.21 

3.84 

10.19 

14.07 

2.68 

0.23 

33.52 

0.00 

100.00 

MEAN 

2.68 

6.86 

13.38 

12.13 

4.48 

8.11 

13.14 

2.20 

0.19 

36.84 

0.00 

100.00 

BOD 

1984 

28.98 

10.76 

15.12 

9.50 

8.15 

11.91 

7.54 

2.27 

0.00 

5.78 

0.00 

100.00 

1985 

28.93 

11.93 

17.24 

7.47 

6.58 

11.41 

8.75 

2.41 

0.00 

5.28 

0.00 

100.00 

1986 

36.56 

5.30 

7.03 

8.81 

4.56 

5.33 

12.12 

6.41 

0.00 

13.88 

0.00 

100.00 

1987 

23.70 

14.10 

18.68 

7.55 

6.99 

11.35 

8.22 

2.55 

0.00 

6.89 

0.00 

100.00 

1988 

31.72 

7.61 

12.92 

6.20 

3.41 

4.01 

11.92 

5.73 

0.00 

16.48 

0.00 

100.00 

1989 

30.05 

9.27 

13.16 

10.34 

6.57 

8.49 

10.19 

3.37 

0.00 

8.56 

0.00 

100.00 

1990 

26.37 

12.21 

17.04 

7.87 

6.89 

11.25 

8.15 

2.78 

0.00 

7.44 

0.00 

100.00 

1991 

28.40 

11.04 

13.19 

9.98 

5.70 

7.09 

9.71 

3.92 

0.00 

10.95 

0.00 

100.00 

MEAN 

28.48 

10.97 

15.21 

8.53 

6.55 

9.79 

9.06 

3.21 

0.00 

8.21 

0.00 

100.00 

SED 

1984 

16.49 

12.77 

21.34 

22.93 

7.16 

19.31 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1985 

14.46 

12.94 

22.36 

21.29 

6.82 

22.12 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1986 

16.14 

11.71 

16.86 

35.96 

5.70 

13.63 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1987 

16.49 

14.87 

23.30 

19.40 

6.59 

19.35 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1988 

21.76 

16.38 

24.37 

24.11 

3.96 

9.42 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1989 

17.53 

11.31 

16.84 

31.38 

6.77 

16.20 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1990 

19.59 

12.18 

20.29 

21.01 

6.73 

20.21 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1991 

16.25 

13.72 

20.18 

30.43 

5.75 

13.68 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

MEAN 

16.99 

13.11 

20.98 

24.07 

6.56 

18.31 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 
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Segment  200 


<• 

FOR 

HTC 

---Pervic 

LTC 

>US . 

PAS 

URB 

HAY 

-X - Impervious-- 

ANHL  RES 

-->  Atnx 
Dep 

3S  Point 
Source 

Sept i c 
Load 

Total 

Load 

NH3 

1984 

3.41 

5.09 

6.22 

8.60 

1.18 

3.74 

4.68 

1.91 

0.51 

64.65 

0.00 

100.00 

1985 

2.30 

2.52 

2.40 

5.92 

0.81 

1.75 

4.92 

1.96 

0.47 

76.96 

0.00 

100.00 

1986 

1.65 

4.36 

4.52 

3.65 

0.53 

1.88 

3.20 

1.83 

0.36 

77.98 

0.00 

100.00 

1987 

2.19 

4.07 

4.34 

5.71 

0.80 

1.96 

4.93 

2.05 

0.50 

73.44 

0.00 

100.00 

1988 

1.27 

3.33 

4.35 

3.80 

0.47 

1.39 

3.33 

1.79 

0.39 

79.91 

0.00 

100.00 

1989 

2.60 

5.91 

7.04 

9.97 

1.26 

4.56 

6.54 

3.54 

0.86 

57.70 

0.00 

100.00 

1990 

3.04 

5.78 

7.73 

10.64 

1.57 

3.59 

7.35 

3.10 

0.75 

56.43 

0.00 

100.00 

1991 

2.87 

3.70 

3.74 

7.23 

0.75 

3.39 

5.19 

2.68 

0.61 

69.83 

0.00 

100.00 

MEAN 

2.39 

4.27 

4.95 

6.71 

0.90 

2.67 

4.86 

2.25 

0.53 

70.42 

0.00 

100.00 

N03 

1984 

18.91 

9.04 

10.95 

22.23 

8.14 

11.57 

0.43 

2.08 

0.56 

0.72 

15.37 

100.00 

1985 

13.27 

9.60 

12.31 

21.79 

8.42 

8.53 

0.54 

2.57 

0.68 

1.01 

21.27 

100.00 

1986 

17.64 

8.48 

9.40 

22.36 

6.08 

8.45 

0.36 

2.68 

0.57 

0,99 

23.00 

100.00 

1987 

34.32 

7.65 

8.12 

18.35 

5.20 

7.81 

0.39 

1.96 

0.49 

0.69 

15.02 

100.00 

1988 

8.82 

12.66 

15.44 

21.13 

5.72 

8.94 

0.38 

2.51 

0,59 

0.94 

22.86 

100.00 

1989 

8.76 

13.24 

16.63 

22.13 

6.54 

10.48 

0.42 

2.58 

0.61 

0.71 

17.90 

100.00 

1990 

6.09 

10.99 

16.01 

24.82 

8.17 

10.74 

0.50 

2.57 

0.64 

0.75 

18.75 

100.00 

1991 

9.39 

8.20 

8.47 

25.67 

7.08 

10.98 

0.41 

2.89 

0.69 

0.97 

25.28 

100.00 

MEAN 

15.62 

9.94 

12.14 

22.10 

6.92 

9.70 

0.43 

2.43 

0.60 

0.83 

19.29 

100.00 

ORGN 

1984 

27.58 

7.31 

8.02 

20.31 

5.34 

3.34 

24.59 

2.04 

0.10 

1.38 

0.00 

100.00 

1985 

29.65 

3.85 

4.53 

15.61 

4.45 

1.83 

34.78 

2.97 

0.14 

2.23 

0.00 

100.00 

1986 

28.64 

2.20 

2.08 

17.11 

4.36 

0.54 

36.39 

5.23 

0.20 

3.27 

0.00 

100.00 

1987 

22.89 

3.88 

3.49 

19.44 

5.10 

0.91 

38.40 

3.45 

0.15 

2.31 

0.00 

100.00 

1988 

27.44 

4.03 

2.46 

11.80 

4.06 

0.39 

40.66 

5.18 

0.22 

3.77 

0.00 

100.00 

1989 

22.99 

4.67 

4.22 

20.73 

5.58 

1.11 

34.68 

3.64 

0.15 

2.19 

0.00 

100.00 

1990 

31.62 

5.35 

5.59 

18.65 

4.96 

2.51 

27.29 

2.42 

0.11 

1.50 

0.00 

100.00 

1991 

31.50 

4,65 

5.00 

10.74 

3.27 

1.82 

34,97 

4.87 

0.20 

2.98 

0.00 

100.00 

MEAN 

27.83 

4.89 

4.97 

17.64 

4.81 

1.86 

32.32 

3.32 

0.14 

2.19 

0.00 

100.00 

TN 


1984 

18.56 

7.90 

9.34 

19.30 

6.08 

7.85 

6.72 

2.07 

0.43 

13.56 

8.20 

100.00 

1985 

15,01 

6.67 

8.26 

16.90 

5.79 

5.42 

8.38 

2.57 

0.51 

19.16 

11.34 

100.00 

1986 

15.59 

6.43 

7.00 

16.87 

4.42 

5.49 

6.04 

2.93 

0.46 

21.42 

13.36 

100.00 

1987 

26.31 

6.34 

6.62 

16.33 

4.39 

5.46 

7.05 

2.29 

0.43 

15,35 

9,43 

100.00 

1988 

9,98 

8,91 

10,54 

15.24 

4.12 

5.63 

6.47 

2.79 

0.48 

22.54 

13.29 

100.00 

1989 

11,14 

10.32 

12.53 

20.22 

5.58 

7.56 

7.63 

3.01 

0.55 

10.22 

11.23 

100.00 

1990 

13.59 

8.64 

11.81 

21.16 

6.31 

7.29 

8.13 

2.65 

0.50 

9.48 

10.42 

100.00 

1991 

13.01 

6.70 

6.94 

19.31 

5.15 

7.72 

7.23 

3.32 

0.58 

14.64 

15.39 

100.00 

MEAN 

15.97 

7.76 

9.19 

18.25 

5.30 

6.59 

7.22 

2.64 

0.49 

15.42 

11.17 

100.00 

P04 

1984 

1.67 

8.26 

12.05 

13.20 

5.00 

13.60 

5.88 

1.92 

0.16 

38.26 

0.00 

100.00 

1985 

1.00 

12.66 

13.25 

12.31 

2.92 

5.84 

6.01 

1.99 

0.17 

43.86 

0.00 

100.00 

1986 

0.29 

11.77 

12.07 

9.14 

1.93 

7.85 

4.50 

2.50 

0.17 

49.79 

0.00 

100.00 

1987 

0.43 

12.61 

13.93 

11.34 

2.91 

9.72 

5.69 

1.99 

0.16 

41.19 

0.00 

100.00 

1988 

0.20 

14.25 

19.54 

5.47 

1.53 

3.53 

4.29 

2.10 

0.16 

48.92 

0.00 

100.00 

1989 

0.51 

16.80 

22.78 

9.29 

3.00 

11.34 

4.81 

1.99 

0.15 

29.31 

0.00 

100.00 

1990 

1.87 

13.25 

16.13 

14.90 

3.97 

12.51 

5.59 

1.95 

0.16 

29.67 

0.00 

100.00 

1991 

0.63 

8.43 

8.83 

12.38 

1.55 

12.69 

4.87 

2.63 

0.20 

47.77 

0.00 

100.00 

MEAN 

0.87 

12.38 

15.11 

11.13 

2.98 

9.80 

5.25 

2.10 

0.17 

40.19 

0.00 

100.00 

ORGP 

1984 

5.20 

15.88 

16.67 

29.66 

5.93 

6.98 

9.56 

2.27 

0.53 

7.32 

0.00 

100.00 

1985 

6.51 

10.03 

11.31 

26.75 

5.90 

4.66 

16.16 

3.94 

0.90 

13.86 

0.00 

100.00 

1986 

5.52 

5.67 

5.17 

29.48 

5.77 

1.37 

16.86 

6.92 

1.28 

21.97 

0.00 

100.00 

1987 

4.47 

9.60 

8.27 

32.25 

6.45 

2.17 

17.00 

4.37 

0.93 

14.49 

0.00 

100.00 

1988 

5.25 

10.14 

5.94 

19.07 

5.21 

0.93 

18.24 

6.63 

1.35 

27.24 

0.00 

100.00 

1989 

4.49 

11.24 

9.75 

33.43 

6.87 

2.59 

14.93 

4.48 

0.91 

11.31 

0.00 

100.00 

1990 

7.98 

13.05 

13.03 

30.52 

6.15 

5.96 

11.85 

3.01 

0.66 

7.81 

0.00 

100.00 

1991 

6.65 

11.87 

12.20 

17.61 

4.27 

4.45 

15.99 

6.36 

1.29 

19.30 

0.00 

100.00 

MEAN 

5.83 

11.86 

11.54 

28.49 

5.95 

4.36 

14.00 

4.11 

0.87 

13.01 

0.00 

100.00 
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Segment  200 


<■ 

---Pervic 

>us - 

-X - Impervious — 

->  Atim 

3S  Point 

Sept  1 c 

Total 

FOR 

HTC 

LTC 

PAS 

URB 

HAY 

ANML 

RES 

Dep 

Source 

Load 

Load 

TP 

1984 

3.13 

11.07 

13.26 

19.61 

5.06 

9.68 

14.21 

1.94 

0.31 

21.75 

0.00 

100.00 

1985 

2.71 

10.69 

11.45 

15.94 

3.64 

4.94 

17.59 

2.45 

0.40 

30.19 

0.00 

100.00 

1986 

1.65 

9.27 

9.35 

13.84 

2.80 

5.55 

14.86 

3.49 

0.45 

38.72 

0.00 

100.00 

1987 

1.62 

10.60 

11.00 

16.74 

3.74 

6.58 

17.29 

2.55 

0.38 

29.48 

0.00 

100.00 

1988 

1.41 

12.21 

14.82 

8.36 

2.31 

2.64 

14.71 

3.04 

0.44 

40.09 

0.00 

100.00 

1989 

1.67 

13.86 

17.15 

15.87 

3.94 

7.86 

15.00 

2.59 

0.37 

21.69 

0.00 

100.00 

1990 

4.07 

12.17 

13.71 

19.76 

4.51 

9.05 

15.13 

2.21 

0.34 

19.05 

0.00 

100.00 

1991 

2.32 

8.77 

9.13 

12.94 

2.22 

9.33 

15.44 

3.51 

0.50 

35.85 

0.00 

100.00 

MEAN 

2.45 

11.25 

12.76 

16.04 

3.73 

7.23 

15.51 

2.61 

0.39 

28.03 

0.00 

100.00 

BOD 

1984 

30.04 

18.23 

12.07 

11.78 

4.82 

9.94 

6.41 

1.85 

0.00 

4.88 

0.00 

100.00 

1985 

35.98 

12.70 

9.87 

10.36 

4.37 

8.33 

9.85 

2.91 

0.00 

5.60 

0.00 

100.00 

1986 

41.10 

5.97 

4.43 

11.93 

4.57 

2.77 

11.01 

5.48 

0.00 

12.76 

0.00 

100.00 

1987 

31.32 

11.43 

7.20 

13.24 

5.33 

4.43 

11.58 

3.61 

0.00 

11.86 

0.00 

100.00 

1988 

35.43 

9.86 

5.89 

7.96 

3.90 

2.92 

11.26 

4.97 

0.00 

17.80 

0.00 

100.00 

1989 

28.72 

13.31 

9.32 

13.15 

5.44 

5.74 

9.75 

3.55 

0.00 

11.02 

0.00 

100.00 

1990 

29.07 

17.00 

11.06 

11.25 

4.65 

10.07 

7.38 

2.27 

0.00 

7.26 

0.00 

100.00 

1991 

37.33 

10.49 

7.04 

7.16 

3.07 

5.35 

9.48 

4.57 

0.00 

15.49 

0.00 

100.00 

MEAN 

32.45 

13.77 

9.29 

11.16 

4.63 

7.20 

8.98 

3.19 

0.00 

9.38 

0.00 

100.00 

SED 

1984 

17.55 

20.57 

18.02 

26.24 

5.99 

11.66 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1985 

22.32 

16.08 

15.36 

29.42 

6.67 

10.14 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1986 

4.40 

15.27 

11.74 

52.62 

10.54 

5.43 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1987 

7.68 

19.16 

13.89 

43.25 

9.53 

6.49 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1988 

2.24 

30.03 

15.09 

37.48 

10.25 

4.92 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1989 

7.59 

20.20 

14.90 

41.15 

9.47 

6.71 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1990 

34.11 

15.02 

12.50 

24.06 

5.55 

8.77 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

1991 

14.86 

24.03 

20.72 

23.01 

5.08 

12.30 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 

MEAN 

18.82 

18.84 

15.39 

30.70 

6.98 

9.25 

0.00 

0.00 

0.00 

0.00 

0.00 

100.00 
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